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ABSTRACT 

Context. The presence of two stellar populations in the Milky Way bulge have been reported recently, based on observations of giant 
and dwarf stars in the inner an intermediate bulge. 

Aims. We aim at studying the abundances and kinematics of stars in the outer Galactic bulge, thereby providing additional constraints 
on formation models of the bulge. 

Methods. Spectra of 401 red giant stars in a field at (/, b) = (0°, -10°) were obtained with the FLAMES-GIRAFFE spectrograph at 
the VLT. Stars of luminosities down to below the two bulge red clumps are included in the data set. From these spectra we measure 
general metallicities, abundances of iron and the a-elements, and radial velocities of the stars. The abundances are derived from 
an interpolation and fitting procedure within a grid of COMARCS model atmospheres and spectra. These measurements as well as 
photometric data are compared to simulations with the Besan^on and TRILEGAL models of the Galaxy. 

Results. We confirm the presence of two populations among our sample stars: i) a metal-rich one at [M/H] ~ +0.3, comprising 
about 30% of the sample, with low velocity dispersion and low a-abundance, and ii) a metal-poor population at [M/H] ~ -0.6 
with high velocity dispersion and high a-abundance. The metallicity difference between the two populations, a systematically and 
statistically robust figure, is A[M/H] = 0.87 ± 0.03. The metal-rich population could be connected to the Galactic bar. We identify 
this population as the carrier of the double red clump feature. We do not find a significant difference in metallicity or radial velocity 
between the two red clumps, a small difference in metallicity being probably due to a selection effect and contamination by the 
metal-poor population. The velocity dispersion agrees well with predictions of the Besanfon Galaxy model, but the metallicity of the 
"thick bulge" model component should be shifted to lower metallicity by 0.2 to 0.3 dex to well reproduce the observations. We present 
evidence that the metallicity distribution function depends on the evolutionary state of the sample stars, suggesting that enhanced mass 
loss preferentially removes metal-rich stars. We also confirm the decrease of a-element over-abundance with increasing metallicity. 
Conclusions. Our sample is consistent with the existence of two populations, one being a metal-rich bar, the second one being more 
like a metal-poor classical bulge with larger velocity dispersion. 
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1. Introduction 

The bulge of the Milky Way galaxy is such a complex system 
that its formation and evolution is still poorly understood. There 
exist two main scenarios for the Galactic bulge (GB) forma- 
tion. The first one, called the "classical" scenario, describes the 
bulge formation through initial collapse of gas at early times 
( |Eggen et al.| |1962|l or through hierarchical merging of sub- 
clumps ( [Noguchi 1999 Aguerri et al. 2001) . In the mono- 
lithic collapse case the bulge formed before the disc and the 



star-formation time-scale was very short (~ 0.5 Gyr; Thomas et 
al. 2005 ). The resulting stars are old (> 10 Gyr) and have en- 
hancements of a-elements relative to iron in a large range of Fe 
abundances, which are characteristic of classical bulges. This in- 
dicates a very fast bulge formation, where SNe la did not have 
time to pollute the gas with a-element-free ejecta. In the hier- 
archical merging the bulge also formed before the disc but on 
a longer timescale (of the order of a few Gyr; NoguchT 1999 1 
and therefore this approach predicts lower over-abundances of 
Qf-elements. 



* Based on observations at the Very Large Telescope of the 
European Southern Observatory, Cerro Paranal/Chile under Programme 
083.D-0046(A). Table [3] is only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-H A/[ 



The second scenario, called the "pseudo-bulge" scenario, de- 
scribes the bulge formation by the buckling of the disc, following 
a disc instability (also called a b ar;,Combes et al. 1990; Raha et| 
al 1991; Norma nlFaLl |T996| [Kormendy & Kenni cutt 2004[ 
|Athanassoula) |2005 1. In this case the bulge forms after the disc 
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and on a much longer timescale than a classical bulge. After the 
bar formation the disc is heated in vertical direction ("Combes' 
& Sanders) |1981 1, giving rise to the typical boxy/peanut shape. 
The bulge formed in this way will be a mixture of disc stars and 
stars formed in situ by gas which is likely to be well mixed by 
the action of the bar In this scenario, the enhancement of the a- 
elements in the bulge stars is predicted to be low, similar to that 
of the inner disc stars in the models. 

It should be noted that, in addition to the classical bulge 
and pseudo-bulge classifications, an alternative formation mech- 
anism has been suggested, namely the "clump-origin bulge". 
In this scenario, first proposed by Noguchi ( 1998| l and recently 
elaborated upon by Inoue 8c Saitoh 1201 1 1, stellar clumps spon- 
taneously form in high gas-density discy galaxies, as shown 
by N-body/SPH simulations, and dynamical friction drags these 
clumps to the centre of the galaxy where they aggregate into a 
bulge-like structure. The models of Inoue & Saitoh ( 20TTll repro- 



duce the observed boxy shape, rapid star formation, and vertical 
metallicity gradients observed in the bulge of the Milky Way. 
Clumpy galaxies that could be the analogues of the primordial 
Milky Way and that are consistent with these models have been 



observed in the high-redshift universe ( Elmegreen & Elmegreen 



2005 



2012 1 



Genzel et al. 2011 1, but also locally ( |Elmegreen et al. 



Each of these scenarios is supported by at least some of 
the observational constraints. The question of formation history 
is crucial and necessary to investigate because our Galaxy is 
a benchmark for understanding the formation of disc galaxies. 
Recently, Babusiaux et al. ( 2010| l and Hill et al. (201 1 1 analysed 
large samples of red clump stars in Baade's window and three 

fields close to the minor axis, atb 4°, -6°, and -12°, which 

revealed the presence of two distinct populations: a metal-poor 
component around [Fe/H] ^ -0.3 dex with a broad distribution 
in [Fe/H], and metal-rich component centered around [Fe/H] » 
+0.3 with a small spread in metallicity. In addition, these two 
populations show kinematical differences: the metal-poor com- 
ponent is compatible with an old spheroid, whereas the metal- 
rich component is consistent with a bar population. Therefore, 
two different formation scenarios have been proposed: A rapid 
formation timescale for the metal-poor component, and a forma- 
tion over a longer timescale driven by the evolution of the bar 
(pseudo-bulge) for the metal -rich component. Similarly, | Bensby| 
et al. ( 201 l| l find two populations in their sample of micro-lensed 
bulge dwarf and sub-giant stars, centred at [Fe/H] ^ -0.6 and 
[Fe/H] ^ +03, with a dearth of stars around [Fe/H] = 0.0. 

Other observations favoured one or the other scenario, which 
prevented a clear picture of bulge formation to emerge. For 
example, [Zoccah et al.| P008j), fr om the observed sample re- 
analysed later by [Babusiaux et al. (2010 1, found a clear metal- 
licity gradient in the bulge; this was interpreted as a challenge to 
the scenario in which the bulge would result solely from the ver- 
tical heating of the bar. On the other hand, recent radial velocity 
studies (jHoward et al.| [20081 12009) [Shen et al.| [20T0l [Kunder[ 
[et al.[ |2012| l find that the radial velocity dispersion (ctrv) of 
bright bulge giants is fully consistent with the pseudo-bulge sce- 
nario. In particular, Shen et aL| ( [2010 ) show that any classical 
bulge contribution cannot be larger than ~ 8% of the disc mass. 
Furthermore, recent spectroscopic studies of the abundances of 
the a-elements as a function of [Fe/H] revealed that at least the 
metal-poor bulge stars are chemically similar to stars in the local 
thick disc (Melindez et al.[[2008[[Alves-Brito et al.[[2010[[Ryde 



A comprehensive summary of the current observational pic- 
ture of the bulge structure, formation, and evolution cannot be 
given in this paper. Instead, we refer to the recent review by [Rich 
( [201 I| l. In any case, it is clear that the current observational evi- 
dence is not sufficient to conclusively constrain the structure and 
formation history of the Galactic bulge, and that more observa- 
tions are required. 

In the past, most of the attention has been paid to the inter- 
mediate bulge, e.g. Baade's window toward (/,/?) - (-1-1°, -3.9°), 
while the outer bulge has been somewhat neglected. This may 
explain why an important feature such as the double red clump 
(RC), which becomes apparent only at \b\ > 5°, was detected 
only recently (Nataf et al., 2010; McWilHam & Zoccali 2010]. 
This feature is interpreted as two over-densities of bulge stars 
at different distances from the sun ( [Mc William & Zoccali] [2010) 
Saito et al. 201 Ij ). It is suggested that the over-densities look like 
a three-dimensional X-structure (Mc William & Zo ccali 2010T). 
A first spectroscopic study of stars located in the two RCs was 



performed by de Propris et al. (2011 1. These authors could not 
discern any difference in kinematics and metal abundance be- 
tween the two RCs, which is in agreement with the interpretation 
that the two RCs represent the same parent population at differ- 
ent distances from the sun. However, these studies were based 
on optical spectra with relatively low resolution and signal-to- 
noise ratio (S/N), from which only an abundance index could be 
measured. An investigation of the double RC at higher spectral 
resolution and S/N is desirable to foster these conclusions. 

In this Paper we present medium-resolution optical spectra 
of ~ 400 bulge stars towards a field at (/, b) = (0°, -10°), i.e. ~ 
1 .4 kpc south of the Galactic plane (assuming a bulge distance of 
8.0 kpc) in the outer part of the bulge. The spectra were obtained 
with the FLAMES spectrograph at the VLT in GIRAFFE multi- 
object mode. Our sample includes stars from the tip of the RGB 
to stars less luminous than the expected bump in the red giant 
branch (RGB) luminosity function (RGB bump), hence it also 
includes the two RCs. This is the first study of the two RCs with 
medium-resolution, high-S/N spectra. We measure from these 
spectra radial velocities, general metallicities, and abundances of 
iron and the a-elements, and combine this information to obtain 
a clearer picture of the structure of the outer bulge. 

Furthermore, we extensively compare our results with pre- 
dictions by models of the Galaxy and use these models to in- 
fer selection biases and the contamination by non-bulge stars. 



Recently, Robin et al. (2012 1 presented a model where two pop- 



[et al.| [2010 1 and the inner disc at Galactocentric distances of 4 
to 7 kpc ( [Bensby et al.[ [2010b[ ). 



ulations co-exist in the bulge region: a bar or pseudo-bulge of 
high metallicity and small scale height, and a "thick bulge" 
or classical bulge with a higher scale height, lower metallic- 
ity, and higher velocity dispersion. This "Besangon model of the 
Galaxy" (BGM) explains well the apparent gradient in metallic- 
ity that is observed along the minor axis, by a variable proportion 
of the two populations of different scale height. Here, we use our 
observed sample as a test case for the new version of the BGM. 

The paper is structured in the following way: The sample 
selection and the observations are presented in Sect. |2] the anal- 
ysis of the data with the help of COMARCS atmospheric models 
and spectral synthesis techniques, as well as with the Besangon 
and TRILEGAL models of the Galaxy, is introduced in Sect. [3| 
Sect. |4] presents and discusses our results on the radial veloci- 
ties, metallicities, and a-element abundances; in Sect. [5] a com- 
parison of the Galaxy models with 2MASS photometry is done; 
finally, conclusions are drawn in Sect. |6] 
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2. Sample and observations 

2.1. Target selection, observations, and data reduction 

Initially, the sample of stars was selected for the study of the evo- 
lution of lithium along the bulge RGB (Lebzelter et al. 2012 1. 
For a detailed description of the sample selection and the obser- 
vations we refer to that paper, here we reproduce only the most 
important points. 

The selection of targets was based on data from the 2MASS 



catalogue (Skrutskie et al. 2006 1 in a 25' diameter circle to- 
wards the direction (l,b) - (0°,-10°), which is the centre of 
the Palomar-Groningen field #3 (PG3). The sample selection is 
illustrated in Fig. [T] which shows a colour-magnitude diagram 
of that circular field. The observed targets (black cycles) were 
chosen to fall close to two isochrones from |Girardi et al.| ( |2000| l: 
Z = 0.004 and age 10 x 10*^ years, and Z ^ 0.019 (which is 
Zq on the scale used by Girardi et al. 2000 1 and age 5 x 10^ 



years. These isochrones were used because they cover well the 
RGB of the bulge. They are not chosen as to reproduce any of 
the GB's properties, in particular not its age spread. A distance 
modulus of 14™5 to the GB was adopted. The chosen targets 
were allowed to have a (7 - Ks)[) colour either bracketed by the 
isochrones, or 0.^02 redder or bluer than them. Dereddening of 
the photometry was done as described in [Lebzelter et al.| ( |2012| l. 
We find a mean reddening of O?' 129 in the 7-band and 0^049 
in the Ks -hand. Besides the colour criterion, only stars were 
included that are fainter than Jo - 9^0 to exclude AGB stars 
above the tip of the RGB, and stars brighter than = 14™5 
to include the RGB bump, which is expected from isochrones 
at 13"8 < 7o < 14™1. This region in the CMD will henceforth 
be referred to as "selection region". Furthermore, all targets were 
excluded that had fewer than two quality flags 'A' in the 2MASS 
JHKs photometry, and targets that had another source within 3" 
that was not fainter by at least 2™0 in the /-band than the tar- 
get itself. Applying all those criteria yielded 514 targets for the 
observations. 

The targets were observed with the FLAMES spectrograph, 
mounted to the VLT UT2 at Paranal observatory, Chile, in the 
GIRAFFE configuration. The grating HR15 centred at 665 nm 
was used, which gave a spectral resolution of 17 000 for a wave- 
length coverage from 644 to 682 nm. The observations were ob- 
tained in service mode in June 2009. Of the initially proposed 
sample, 401 spectra had sufliciently high quality to measure their 
radial velocity (RV, see below). 

2.2. Magnitude liistograms and tiie two RCs 

Fig. |2] shows histograms of the 2MASS magnitudes of all stars 
in the selection region (527 stars, dotted line), and of the stars 
for which the RV could be determined (401 stars, solid line). 
At intermediate magnitudes (13™0 < 7o < 11™5), spectra of 
only about 40% of the selected targets are available because one 
fibre configuration for the stars in this brightness range was not 
observed. The reduced fraction of actually observed stars will be 
taken into account in the analysis (Sect. |4TTTT] i. 

Also discernible in Figs. \l\ and |2| are the two red clumps 
(RCs) rece n tly ide ntified by Nataf et al. (2010) and Mc Wilhaml 
|& Zoccali| ( |2010| l. These structures were not known to us in 
the design phase of the programme. Initially, the fainter RC at 
Jo ~ 13™9 was interpreted as to be the RGB bump, and only 
the brighter one at Jo ~ 13™2 as the red clump. Inspecting 
Fig. |2] we defined membership of stars in the two RCs as fol- 
lows. Bright RC: 13™65 > 7 > 13'."05, 13^05 > H > 12'"40, 



Z=0.004, 10x10'' yrs 
Z=0.019, 5x10^ yrs 




1.4 



Fig. 1. De-reddened 2MASS colour-magnitude diagram of the 
25' diameter field around {l,b) - (0°,-10°). Targets observed 
in the present programme are plotted as black circles, those that 
have not been observed as black dots. The two RGB isochrones 
from Girardi et al.| ( |2000) used for the target selection, with ages 
and metallicities as indicated in the legend, are shown as dashed 
and solid lines. The isochrones are truncated at the tip of the 
RGB. The two RCs are the over-densities at Jo ~ 13™2 and Jo ~ 
13?9. 
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Fig. 2. Histograms of the observed (not de-reddened) 2MASS 
magnitudes of our sample stars. The dotted line shows the his- 
tograms of all stars in the selection region (527 stars, including 
those that did not fulfil the quality flag criteria to be selected 
as target), the solid line shows the histogram of the 401 stars 
for which the radial velocity could be determined. The vertical 
dashed lines indicate the magnitude limits that were adopted to 
define the memberships in the two RCs (Sect. 2.2 1. The double 
RC is clearly visible in all three bands. 



and 12'7'95 > K > 12™30; faint RC: 14™50 > 7 > 13'?^65, 
13™80 > H > 13'J'05, and 13'?70 > > 12?95. We found 96 
and 150 stars, respectively, which fulfill these criteria. Because 
it is impossible to ascribe with certainty a given star to either 
the front or the back arm of the X-structure, the uncertainty of a 
bulge star's absolute magnitude is about 0'."35. 
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Table 1. Parameters of the COMARCS model grid. 
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Fig. 3. Histogram of the total proper motion of our sample stars. 
Stars with a total proper motion larger than 20 mas/year are de- 
fined as foreground candidates. Note that four sample stars with 
a total proper motion larger than 60 mas/year fall outside the plot 
range. 
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ground candidates. Probably there are still many genuine bulge 
stars among these foreground candidates, which is why we re- 
frain from excluding these stars from the analysis. 



3. Analysis 

The spectra were analysed with the help of model atmospheres 
and synthetic spectra, while the results are interpreted with mod- 
els of the MiUcy Way galaxy. Both these types of models and the 
associated analysis are described in the following. 



2.3. Foreground contamination from proper motions 

An important issue for this kind of study is the foreground con- 
tamination of the sample. The proper motion can be used to 
identify possible foreground stars. We therefore searched the 
recently published Southern Proper Motion Program catalogue 
version 4 (SPM4; ^Girard et al. 201 1 1 for sample stars with large 
proper motion. All but ten of our sample stars were found in the 
SPM4 catalogue. In principle, foreground stars could be iden- 
tified as those stars whose total space velocity clearly exceeds 
the maximum measured RV, assuming that they are placed at 
the distance of the GB (~ 8 kpc). The maximum RV found in 
our sample is about 300 km s"', which corresponds to a (tangen- 
tial) proper motion of 7.91 milli arc-seconds (mas) per year at 
a distance of 8 kpc. Unfortunately, the mean combined error on 
the proper motion quoted in SPM4 (6.3 mas/yr) is only barely 
smaller than this value, so with this catalogue it is impossible to 
probe in detail the proper motions out to the GB. Applying the 
criterion described above necessarily would return a mixture of 
bulge and disc foreground stars. Figure [3] shows the distribution 
of the total (combined) proper motions of our sample. The distri- 
bution essentially consists of a broad peak centred on ~ 7 mas/yr 
and a long tail of higher proper motion stars. The tail seems to 
start at ~ 20 mas/yr, which is why we decided to assign those 
stars to the foreground whose combined proper motion is larger 
than this value. This is also approximately the number that one 
gets when two times the mean combined error (2cr) is added to 
the most common total proper motion. Applying this criterion, 
we found 50 (12.5%) foreground star candidates, which roughly 
agrees with the fraction of non-bulge stars estimated from the 
Galaxy models (see Sect. 3.2.3 1. The foreground star candidates 
are marked with an asterisk in Table [3] One star in our sam- 
ple (#300) has a particularly large proper motion of more than 
200 mas/yr, hence it is probably a nearby K dwarf. In addition, 
#143 is probably a foreground M dwarf. 

Unfortunately, there is reason to believe that there are large 
(systematic) errors in the SPM4 catalogue because the dou- 
ble RC as well as the two peaks in the metallicity distribution 
(Sect. |4.1.2| i, two features which are believed to be related to the 
GB rather than to the disc, are also found among these 50 fore- 



3.1. COMARCS atmosptiere models and COMA syntlietic 
spectra 

For the analysis of the spectra we used model atmospheres calcu- 
lated with the COMARCS code (Aringer et al.!l 2009| l, and spec- 
tra based on these models synthesised with COMA (Copenhagen 
Opacities for Model Atmospheres; [Aringer et al. 2009 | l. 

A determination of the stellar temperature and surface grav- 
ity with a classical equivalent width analysis using only the 
FLAMES spectra was not possible due to the relatively small 
spectral range and the dominance of TiO lines in the cooler tar- 
gets. Instead, the temperatures of the stars were determined from 
their 2MASS {J - Ks )o colour and the effective temperature cal- 



ibration based on COMARCS models presented in Lebzelter et 
al. (2012). Hence, the calculations to derive temperatures are 
consistent with the spectral synthesis applied in this work to 
obtain abundances in the stars. The surface gravities (g) of the 
stars were obtained from the Jq magnitudes of the stars and 
the isochrones in Fig. [T] assuming that the stars belong to the 
bulge RGB. The uncertainty in the distance modulus of ±0^35 
(Sect. 2.2 1 introduces an uncertainty in log g of +0. 16 dex, which 
however has only a minor impact on the uncertainty in the abun- 
dances. Furthermore, a generic micro-turbulence of 2.5 km s"', 
typical of red giant stars, a C/O ratio of 0.3, and a generic over- 
abundance of the ff-elements (O, Ne, Mg, Si, S, Ar, Ca, and Ti) 
of -H0.2dex were adopted for the model calculations. This over- 
abundance was adopted because previous spectral abundance 
studies of bulge stars showed a general enhancement of the 
Qf-elements of this order of magnitude (e.g. |Alves-Brito et al.| 
|20T0| . 

With the temperature, \ogg, and micro-turbulence fixed in 
this way, a grid of model atmospheres and associated synthetic 
spectra was calculated. To keep the grid within reasonable size 
limits, the RGB was sampled by five different values of log g. For 
each of these bins, models at various temperatures, in a range ap- 
plicable for this log g bin, were calculated. For each combination 
of log g and T^ff, model atmospheres at five different metallici- 
ties [M/H] = 1.5, -1.0, -0.5, 0.0, and -1-0.5 were calculated. The 
parameters of the model grid are summarised in Table [T] 

In the next step, spectra were synthesised for all the models 
in the grid and convolved with a Gaussian to the resolving power 



4 



Uttenthaler et al.: The structure and formation of the Galactic bulge 



of the observed spectra (R - A/AA - 17 000). The most impor- 
tant hne lists used are the one for TiO ( Schwenke! 1998 1, CN 
(J0rgensen 1997), and the atomic lines (VALD; Kupka et al. 



999| l. This last list was checked and improved by comparing 
a model spectrum of Arcturus (adopting the stellar parameters 
and abundances found by |Ryde et al. 20101 with the observed 
high-resolution Arcturus spectrum (Hinkle et al.[|200 0). The so- 
lar abundances listed in Caffau et al. (2008| were adopted as ref- 
erence scale. The solar metalhcity on that scale is Zq - 0.0156. 

We emphasise that a particular strength of our method to 
measure abundances is the consistency of the temperature deter- 
mination as well as model atmosphere and synthetic spectra cal- 
culation, for which the same opacity tables and radiative transfer 
routines were applied. 



3.1.1. Radial velocity measurement 

The first quantity that was measured from the spectra was the 
heliocentric radial velocity. This was done with the help of a 
cross-correlation technique over the whole available wavelength 
range, using a synthetic spectrum as template. At this step, a 
model spectrum with similar temperature as the star under in- 
vestigation was used for the cross-correlation. The abundance of 
the model spectrum was not adjusted for each star, we simply 
adopted solar abundances. The RV was successfully determined 
for 401 sample stars. We estimated a typical error of ~ 0.5 kms"' 
on the measured RV from a Gaussian fit to the peak in the cross- 
correlation function. 
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Fig. 4. Typical examples of the fit qualities achieved by our grid 
interpolation. Observed spectra are plotted as solid black lines, 
best-fit interpolated model spectra as red lines. From bottom 
to top: star #172 {T^ff = 4444 K, log^ = 2.1, [M/H]=-i-0.06, 
Xmin = 0.0446), star #355 (T^ff = 4547 K, log g = 2.5, [M/H]=- 
OM,Xmm = 0.0303), and star #191 (r^ff = 4569 K, log g = 2.1, 
[M/H]=-1.05,A'mm = 0.0166). The specti-a of #355 and #191 ai-e 
shifted vertically by 0.25 and 0.5 for clarity. A few prominent 
metal lines are identified. 



3.1 .2. Determination of the metal abundance 

To determine the general metal abundance from the spectra, we 
employed a simple interpolation and fitting procedure. First, the 
value for \ogg was chosen according to the Jq magnitude of 
a given star (see Table [TJ- Then, all model spectra with that 
log g and a given metallicity were interpolated to the tempera- 
ture of the star using spline polynomials. This interpolation was 
done for every point of the wavelength vector For the subse- 
quent fitting procedure, the observed spectrum had to be nor- 
mahsed. For the hotter stars (reff > 3980 K, corresponding to 
(J - Ks)o < 0.93) this was done by adjusting the observed flux 
to the model spectrum flux around a few continuum points. In 
the spectra of the cooler stars, lines of TiO are appearing or even 
dominating the spectrum, which prohibit a safe definition of the 
continuum; therefore, for these stars the observed spectrum was 
normalised as to have the same median flux as the model spec- 
trum over the used wavelength range. The wavelength range was 
limited to 649 - 680 nm (vacuum wavelength). Finally, the metal 
abundance was determined by a minimisation method. The 
figure of merit to minimise was chosen to be 



N 



;^([M/H]) 
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where is the number of wavelength points, and /obs,, and 
/model,/ are the observed and synthetic model flux at wavelength 
point /. Hence, x is the mean difference between observed and 
synthetic model flux, in units of the observed flux. The IDL rou- 
tine amoeba. pro was used to find the metal abundance [M/H] 
where x reaches a minimum. The general metal abundance was 
successfully determined for 383 sample stars. 

A few examples of the typically achieved fit qualities are 
displayed in Fig. |4] This figure shows a zoom-in on the group 



of metal lines located between 649 and 65 1 nm. An advantage 
of our fitting method is that it also works well at relatively low 
S/N because all the information in the whole wavelength range 
is taken into account. A good example for this is star #172 in 

Fig.g 

An illustrative check of our metallicity determination is dis- 
played in Fig.|5] which shows a CMD of our sample stars where 
the colour of the plotting symbol codes the metallicity of each 
star Metal-rich stars tend to be more common at the red edge 
of the selection region, whereas the metal-poor stars tend to be 
more abundant on the blue edge. This is what can be expected 
from the evolution along the RGB: at a given brightness, a metal- 
rich star will be cooler, hence redder in (J - Ks){), than a metal- 
poor star. This suggests that our methods for temperature and 
metallicity determination work well. Additional sub-structure 
around the double RC is hidden by the high density of stars in 
that figure; it will be discussed in more detail in Sect. 4.3 



There are a few sample stars for which the fitting rou- 
tine extrapolated beyond the abundance range covered by the 
COMARCS model grid. The ten metal-poor outliers clearly 
have very weak metal lines, indicating a very low metallicity 
([M/H] < -1.5). With such a low metallicity, these stars pos- 
sibly belong to the halo population. The two metal-rich outliers 
are found at the bright and cool end of the sample. 



3.1 .3. Determination of the Fe and a-abundance 

An estimate of the abundances of iron and the a-elements can 
be achieved with a procedure very similar to the one above. The 
only difference here is the non-uniform weighting in wavelength. 
Weight is put only to those wavelength points that are influenced 
by line absorption by an or-element or Fe, respectively. Hence, 
the X value to minimise was altered to 
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Fig. 5. Colour-magnitude diagram of our sample stars, with the 
metallicity colour-coded as indicated in the legend. The sample 
was divided at [M/H] - -0.4 into two approximately equally 
large parts. In particular on the upper RGB, metal-rich stars are 
preferentially located at the red edge of the selection region, 
whereas metal-poor stars are preferentially at the blue edge. 
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where W, is the weight function. It is a binary function with 
value 1 for wavelength points with absorption by Fe or an a- 
element of more than 2% of the continuum flux, and oth- 
erwise. The weight function was calculated from a synthetic 
spectrum based on the model atmosphere with Tefr = 4450 K, 
loggr - 2.1, and [M/H] = -0.50, taking into account only the 
a-elements (O, Ne, Mg, Si, S, Ar, Ca, and Ti) or Fe, respec- 
tively. The different a-elements will not contribute equally to 
this mean abundance. In particular, O, Ne, and Ar have no note- 
worthy lines in our wavelength range. For the other elements, 
the relative weight on the average c-abundance is: Ca : Ti : 
Si : Mg : S = 1.000 : 0.751 : 0.400 : 0.027 : 0.009. 
Hence, the a-abundance is dominated by the individual abun- 
dances of Ca, Ti, and Si. The interpolation was done between 
the grid points at [a/R] = (-1.3, -0.8, -0.3, -hO.2, -hO.7) and 
[Fe/H] = (-1.5,-1.0,-0.5,0.0,-^0.5), respectively. The fitting 
routine did not converge for all stars. Eventually, a-abundances 
were determined for 368 stars, and Fe abundances for 363 stars. 
The determination of the iron and a-abundance failed in partic- 
ular for many of the cool, TiO-dominated stars. The Fe abun- 



dance is a very good tracer of the general metal abundance: the 
mean of [M/H]-[Fe/H] is -0. 13dex, with a standard deviation 
of 0.12 dex; the standard deviation decreases to only 0.07 dex 
among the sample stars fainter than Jq = 12™7. 

All measured abundances and radial velocities can be found 
in Table [3] in the on-line Appendix. 

3.1 .4. Systematic errors 

A check of our relatively simple methods for the metallicity and 
abundance determination was done by applying it to the spec- 
trum of Arcturus. For this star, our method yields a metallicity 
of [M/H] = -0.72, and [Fe/H] = -0. 92. If dedic a ted m odel 
atmospheres with the parameters from |Ryde et al.| ( 2010) 1 are 
used (Teff - 4280 K, \ogg - 1.7), i.e. if the grid interpolation 
in temperature is avoided, [M/H] increases to -0.69. Models 
without an a-enhancement (i.e. scaled solar abundances) yield 
[M/H] = -0.64. A search on the Simbad databas^ shows that 

most studies find [Fe/H 0.50 for Arcturus. 

To back-up the existence of a systematic abundance offset. 



we made a blind test with four stars from the sample of Gonzalez 
| etal.| ( |2011a[ ), which cover a range of almost 1 dex in [Fe/H] . The 
spectra of the stars in the FLAMES HR15 setting were kindly 
provided by O. Gonzalez (private communication) and their stel- 
lar parameters as determined in |Zoccali et al.| p008) were com- 
municated, but not their individual abundances. The same fitting 
routine as for our sample stars were appUed to the spectra. The 
results were then compared to those of |Gonzalez et a l.'(201 la|l. 
The Fe abundances derived with our method were systematically 
lower by 0.19 to 0.29 dex than those of Gonzalez et al. ( |2011a i. 

Systematic abundance offsets of 0.2 to 0.3 dex between 
equivalent width measurements and;\f^ minimisation techniques 
(in the sense that the minimisation gives lower abundances) 
have also been reported by the RAVE survey team (Zwitter et ak] 
[2008 Boeche et al. 201 1 1, in particular for giant stars. No expla- 
nation as for the possible origin of this offset is given in the liter- 
ature, neither did we attempt to track down its cause. One expla- 
nation could be that the classical equivalent width measurement 
uses only relatively few lines, selected to be unblended, whereas 
the minimisation uses many more lines, even blended and 
partially saturated ones. While the reasons are only speculative 
at this point, we decided to add 0.2 dex to all abundances deter- 
mined here. This does not only yield better agreement with the 
comparison stars of Gonzalez et al. (201 la) and Arcturus, but, as 



we will see below, gives very satisfactory agreement with other 
studies of abundances in GB stars. The thereby corrected abun- 
dances will be used henceforth in the analysis of our data set 
and are also reported in Table [3] Despite the systematics in the 
abundance scale, we are nevertheless confident that our abun- 
dance determination is reliable in the relative sense within the 
sample. Furthermore, we decided to base our analysis mainly on 
the [M/H] measurement because it is available for a few more 
stars, in particular more of the bright stars. The results would 
remain qualitatively unchanged if [Fe/H] was used instead. 

3.1.5. Uncertainty estimate 

The uncertainty in the determination of the metallicity was es- 
timated for six representative sample stars covering a range in 
temperature, metallicity, brightness, and S/N ratio. This was 
done mainly by varying the free parameters within their uncer- 
tainties, one at a time while keeping the others constant, and then 



http://simbad.u-strasbg.fr 
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repeating the fitting process. Uncertainties arising from the de- 
termination of temperature and log g, as well as in the continuum 
placement, were taken into account. 

The uncertainty in the effective temperature was determined 
from the uncertainty on (J - Ks )o by adopting the uncertainties 
in the J- and Ks -hand magnitudes given in the 2MASS cata- 
logue. For the bright (cool) stars this uncertainty is between 60 
and lOOK, while for the faint (hot) ones it is between 130 and 
200 K. However, this does not mean that the temperature has a 
lower impact on the [M/H] determination for the bright and cool 
stars, just on the contrary: because their spectra are dominated by 
the very temperature-sensitive TiO bands, a change in tempera- 
ture has the largest impact on their metallicity determination. 
Uncertainties up to 0.27 dex from the temperature uncertainty 
alone were found, but for most stars it is in the range 0.09 to 
0.18 dex. 

The uncertainty coming from the log g estimate was deter- 
mined by doing the interpolation and fitting with a group of 
model spectra adjacent in \ogg to the one that would be at- 
tributed to a given star based on its magnitude. We found that 
the \ogg estimate contributes 0.06 - 0.10 dex of uncertainty in 
the [M/H] determination. This should already include the uncer- 
tainty coming from the fact that we interpolate model spectra 
to the temperature of the stars, rather than using model spectra 
calculated for the exact temperature of every star. 

Finally, the uncertainty introduced by the continuum place- 
ment was estimated by varying the normalisation factor within 
its uncertainty range, derived from the scatter of the pixels in the 
continuum points. This uncertainty mainly depends on the S/N 
of the spectrum, and contributes another 0.02 - 0.04 dex to the 
uncertainty in [M/H]. Thus, for most stars the combined uncer- 
tainty in [M/H] is in the range 0.1 1 - 0.21 dex, but it can be as 
high as 0.30 dex for the cool stars dominated by TiO absorption. 
The uncertainties in the iron and a-element abundances are very 
similar. 



3.2. The Besangon and TRILEGAL models of the Galaxy 
3.2.1. The Besangon model 

The Besangon Galaxy Model (BGM) is based on assumptions on 
the scenario of formation and evolution of four main stellar pop- 
ulations of the Milky Way ( Robin et al.||2003) . It allows to sim- 
ulate the stellar content in any given line of sight, and for each 
simulated star the photometry, kinematics, and metallicity are 
computed. For each population, the main assumtions are: a star 
formation-rate history, age, and initial mass function, which al- 
low to generate a distribution function in absolute magnitude, ef- 
fective temperature, and age of the stars. These assumptions are 
confined by observations, such as the stellar content in the solar 
neighbourhood, remote star counts, and photometry in the visi- 
ble and near-infrared. Density laws are assumed and controlled 
by dynamical principles ( jBienayme et al.[ [1987 ) for each popu- 
lation and are tested by means of photometric star counts. The 
model has been extensively compared with 2MASS data which 
allow to constrain the thin disc parameters quite well, such as 
the scale length, scale height, the inner disc hole, and the warp 
(Reyle et al. 2009). However, some model parameters are not 
well constrained yet, such as for the thick disc population. When 
a parameter is not yet well constrained from available data, such 
as metallicity gradients for example, a mean value from the lit- 
erature is assumed, or a value expected from standard scenario 
of formation. For example, thick disc parameters are still con- 
troversial: its density parameters such as the local density and 



the scale height are degenerated when obtained by model fitting 
from high galactic latitude fields ( Robin et al.[ 2012 1. However, 
in the direction considered here, the thick disc population re- 
mains marginal. We conservatively use here the thick disc pa- 
rameters determined from high latitude counts by Reyle & Robin 
( |200l] l. 

Once density functions and distribution functions are ob- 
tained, the simulated stars are generated in the whole Galaxy di- 
vided in volume elements. A number of improvements have been 
made to the model since the one presented in Robin et al. ( 2003 i. 
The mean metallicity and dispersion about the mean for each age 
bin in the thin disc are now taken from Haywood (2008), rather 
than fromlTwarog|( 1980|). The model also includes a 3D extinc- 
tion map (Marshall et al. 2006| l and photometric errors are added 
for comparison with observations. 

Recently, a new model of the bulge/bar region has been pro- 
posed by Robin et al. (2012), which we will use here to interpret 
and compare to our observed data. The new model includes a bar 
and a "thick bulge". In the following, this latter population will 
be referred to as the bulge when discussing the BGM. The bar is 
the most massive component, which dominates the stellar con- 
tent at low latitudes. The bulge is longer and thicker and gives a 
contribution at intermediate latitudes where the bar starts to be 
less prominent. The main bar population has a mean metallicity 
close to solar, while the bulge population is found to be more 
metal-poor, with a metallicity of the order of -0.4 to -0.5 dex. 
The mean metallicties and dispersions of these two populations 
have been roughly estimated from the observed MDF by Zoccali 
et al. ( 2008 1. It should be better constrained in the future using 
larger samples and more lines of sights in the inner and outer 
bulge. 

A large hole in the central region causes the thin disc to reach 
its maximum density at about 2.5 kpc from the Galactic centre. 
The kinematics of the thin disc follow empirical estimates from 
Hipparcos (Gomez et al. 1997), while for the thick disc it is 
taken from proper motion data from ,Ojha et al.|(|1996), and for 
the halo from |Norris et al.| ( fT985| ), as in [Robin et al.| ( |2003| ). hi 
the version presented here, the bar kinematics are taken from the 
dynamical model of |Fux]([T999]l. The bulge is found to have a 
larger velocity dispersion, close to the one of the halo (Robin et 
al., in preparation). 

The fact that the bar population has a smaller velocity dis- 
persion, a higher metallicity, and a smaller scale height than the 
bulge allows to explain well the vertical gradient of metallicity 
seen in bulge fields along the minor axis (Zoccali et al. 2008| l. 
|Babusiaux et al.| ( [2010| l already suggested that the origin of the 
metallicity gradient could be a varying mix of these two popu- 
lations. The bar and bulge populations in the BGM do not have 
an intrinsic metallicity gradient. This model also explains well 
the presence of double clumps at medium latitudes, as seen in 
Nataf et al.| ( |20T0l i, |Mc William & ZoccaE] ( |20T0l l, and |Saito et| 



al. (201 1 1. This derives from the fact that the bar flares from in- 



side to outside by about 30%. This flare can be due to resonances 
in the bar that trap thin disc stars and scatter them to higher scale 
height. 

3.2.2. The TRILEGAL model 



The TRILEGAL model was developed by |Girardi et~aL] ( |2005) 
and is a population synthesis code for simulating the stellar pho- 
tometry of any field in the Milky Way galaxy. The model has 
proven to well reproduce number counts, amongst others, in all 
three pass-bands of the 2MASS catalogue, with errors smaller 
than ~ 30 per cent ( |Girardi et al.| |2005| ). The bulge component 
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was introduced in TRILEGAL only later by |Vanhollebeke et al.| 
( [200 9'), who also calibrated its stellar population and metallicity 
distribution using photometry from the 2MASS and OGLE-II 
surveys of red clump stars. The input MDF for the bulge compo- 
nent in TRILEGAL was obtained by |Vanhollebeke et al.| ( [2059) l 
with the help of a fitting procedure to the RC star distribution, 
by applying shifts to the photometric MDF of |Zoccali et al.| 
(2003 ). The best fit was obtained with a shift of +0.3 dex. Note 
that Vanhollebeke et al. (2009) used fields that are considerably 
closer to the Galactic plane than our field, at most 5.91° south of 
the plane. 

Their best-fit bulge model is provided as the default in the in- 
teractive web interface to TRILEGA1_0 which we also used here 
to simulate the stellar content of our FLAMES field. Other im- 
portant ingredients that we adopted included: a binary fraction of 
0.3 with mass ratios between 0.7 and LO; an exponential extinc- 
tion law with an extinction value at infinity equal to the one used 
to deredden the observed photometry (Cr526 in the V-band); a 
thin and a thick disc, both with squared hyperbolic secant den- 
sity distributions; and a halo with an oblate r'^"* spheroid density 
distribution. For all these components, the default parameters of- 
fered on the web interface were adopted. The TRILEGAL model 
does not provide for the kinematics of its simulated stars, hence 
a comparison with our observations was not possible in this re- 
spect. 

A circular field of 25' diameter centred on the point (/, b) - 
(0°,-10°), equal to the area from which the observed targets 
were selected, was simulated with both models. An exception 
here is that for the BGM we also calculated ten simulations 
with identical parameters, to decrease the Poisson noise mainly 
for the comparison with the observed RVs. For all subsequent 
comparisons between spectroscopic observations and the Galaxy 
models, the same photometric selection criteria were applied to 
the models as to the observed sample (Sect. 2. 1 1. In Sections |4] 
and|5] we will confront the model predictions with our observa- 
tions in terms of kinematics, abundances, and photometric prop- 
erties. 



3.2.3. Inferring metallicity biases, foreground contamination, 
and mass distribution with the Galaxy models 

The Galaxy models can be very useful to infer a number of im- 
portant properties of the observed sample, in particular metalli- 
city biases imposed by our photometric selection criteria, (fore- 
ground) contamination by non-bulge stars, and the distribution 
of initial masses. 

The location of a bulge star in a CMD, and hence its pro- 
bability to fall in our selection region, is determined by a com- 
plex interplay of metallicity, distance, age, reddening, and pho- 
tometric errors. The ratio of the input MDF for the bulge (and 
bar) component to the MDF of the stars selected with our pro- 
cedure reveals any such metallicity biases. An inspection of this 
ratio for the BGM simulation shows that there exists a small bias 
against intermediate metallicities ([M/H] 0.4). This is prob- 
ably caused by the different distance distributions of the bulge 
and the bar component in the BGM: the bulge component, which 
contributes most stars in this intermediate to low metallicity, is 
spread out between 2 and 14kpc distance from the sun, whereas 
the bar, which is centred around solar metallicity, is much more 
concentrated around 7 kpc distance. Many of the nearby bulge 
stars thus do not fall in the selection region, which causes this 
small bias against intermediate metallicity stars. In fact, an in- 



spection of the distance distribution shows that our selection 
procedure is relatively efficient in picking up stars at distances 
between 5 and 1 1 kpc. 

The TRILEGAL model includes a single bulge component, 
with a distance distribution centred on 7 kpc. Hence, the metal- 
licity bias inferred from this model somewhat differs from the 
one found from the Besangon model. The ratio of input MDF 
to the MDF in the selection region is relatively flat in a broad 
range from [M/H] ~ -0.4 to -1-0.5, but gradually decreases to 
lower metallicities. This is also what one would expect from 
a bulge population with a relatively small distance spread: be- 
cause low-metallicity stars evolve along the RGB at relatively 
blue J - Ks colours, in particular on the upper RGB, they would 
be missed by our selection procedure (along relatively metal- 
rich isochrones) with a higher probability. Consistent with this 
expectation, we find in both models that the metal-poor stars in 
the selection region are preferentially at larger distances than the 
metal-rich ones because the blue, metal-poor stars are shifted 
into the selection region when they are at large distances. 

Both the Besangon and TRILEGAL model allow to assign 
a simulated star to one of the stellar components of the Milky 
Way, by which means we estimated the fraction of genuine bulge 
stars in our sample. The selection of stars from the Besangon 
model yielded 688 stars, of which 333 (48.4%) belong to the 
thick bulge, 216 (31.4%) to the bar (hence 79.8% of bulge/bar 
stars in total), 72 (10.5%) to the thin disc, 64 (9.3%) to the 
thick disc, and 3 (0.4%) to the halo. On the other hand, in the 
TRILEGAL model 509 stars fall in the selection region, of which 
479 (94.1%) belong to the bulge, 21 (4.1%) to the thin disc, 9 
(1.8%) to the halo, and none to the thick disc. From this esti- 
mate we conclude that a large fraction (> 80%) of our sample 
stars are genuine bulge stars. Most of the foreground contami- 
nation comes from the thin disc, while the halo contamination 
should be very small. The fraction of non-bulge stars agrees 
reasonably well with the fraction of foreground star candidates 



selected on basis of their proper motion (12.5%) in Sect. 2.3 



http://stev.oapd.inaf.it/trilegal 



Most notably, applying the same proper motion selection crite- 
rion of > 20 mas/yr to the simulated stars in the Besan§on model 
retrieves 52 foreground star candidates, in excellent agreement 
with the number we found from the SPM4 catalogue. Of these 
52 stars, 5 1 are thin disc stars and one is a thick disc star. 

The distribution of initial masses of the simulated stars 
turns out to be relatively sharply peaked in both models. In the 
Besangon model the mean mass is 1.O5M0 with a standard devi- 
ation of 0. 1 IMq, whereas in TRILEGAL the mean is at 1 . 13M0 
with a standard deviation of O.14M0. Genuine bulge stars have 
an even narrower distribution: 1 .07 + O.O5M0 in the Besangon 
model and 1.15 + O.O9M0 in the TRILEGAL model. The most 
massive stars selected from the simulated data have 1.39 and 
1.59M0 in the Besangon and TRILEGAL model, respectively. 
Hence, we are confident that our sample contains mostly low- 
mass (M ~ I.IM0) red giant stars located in the GB. 



4. Results and Discussion 

4.1. The observed metallicity distribution function 
4.1 .1 . Whole sample 

The MDF of our sample stars is displayed in Fig. |6] Two ver- 
sions are plotted: The raw MDF (solid line), and a MDF cor- 
rected for the incompleteness in some brightness (and hence 
colour) ranges in our sample. To correct for this, we appUed 
a bootstrapping method similar to method ii) used by |ZoccaIi| 
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Fig. 6. Metallicity distribution function of all 383 sample stars 
with determined metallicities. The solid line is the raw MDF, 
while the dotted line is an MDF corrected for the non-uniform 
fraction of observed vs. available targets as a function of magni- 
tude, using a bootstrapping method. Poissonian noise error bars 
are also shown for the raw MDF. 



et al. (2008 1. In this method, stars are drawn randomly until 
a flat ratio of observed to available stars is reached. IZoccali etl 
^ ([2008") considered this ratio as a function of colour. Instead, 
we chose to consider it as a function of brightness because our 
targets were grouped for observations according to their bright- 
ness, where one of the scheduled settings for the intermediate 
brightness group was not observed. However, since our selec- 
tion region diagonally crosses the CMD, any under-sampling 
in colour bins will be effectively corrected for by this method. 
Two-hundred such randomly drawn MDFs were created and av- 
eraged, the result of this procedure is plotted as dotted line in 
Fig. |6] The corrected MDF is only marginally different from the 
raw MDF, hence we conclude that our sample well represents 
the metallicities of stars in the whole selection region. 

The MDF is broad with a clear peak at [M/H] ~ -0.55. 
The median, mearj^ and standard deviation of the MDF are 
-0.41, -0.34, and 0.51 dex. This mean metallicity agrees well 
with what |Zoccali et a l. (2008) found: for fields at b - -6° 
and ^ = -12° they found <[Fe/H]) = -0.21 and -0.29, respec- 
tively. Interpolating between these fields to b = -10° yields a 
mean of -0.26. Note, however, that |Zoccah et al.| ('2008) esti- 
mated a large contamination by thin and thick disc stars in their 
sample at b - -12°, which may have a considerable impact on 



the mean metallicity of this field. Whereas Zoccali et al. 



( [2008) 1 



used a box in an / vs. V - / colour-magnitude diagram to select 
their targets, we used a selection region along isochrones, a strat- 
egy that is probably much more discriminant against thin and 
thick disc contaminants. Johnson C. et al. ( 201 1\ derived spec- 
troscopic iron abundances of RGB and RC stars in two fields 
at (l,b) = (0°,-8°) and (l,b) = (-l°,-8.5°). The mean [Fe/H] 
found for these fields are -0.34 and -0.27, respectively. These 
authors also derive photometric metallicities for the latter field 
and find a median of [Fe/H] - -0.34. Also [Gonzalez et al.| 



(201 lb I derived photometric metallicities along the GB minor 
axis and find <[Fe/H]> = -0.36 in a 30' x 30' field around 



^ Doing the averaging on a linear scale, which is the mathematically 
correct way, the mean is [M/H] = -0.05, converted back to the loga- 
rithmic scale. 
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Fig. 7. MetalHcity distributions of bright (/q < 12"7, dashed 
line) and faint (Jq > 12™7, solid line) sample stars. The metal- 
poor and metal-rich populations are more clearly separated 
among the faint sample stars. 



b — -8°. These values all agree nicely with one another, and 
with our result for the field atb - -10°. 

The MDF in Fig. [6] exhibits a relatively sharp cut-off at 

[M/H] 1.0. This could be related to the use of an isochrone 

with Z = 0.004 (corresponding to [M/H] = -0.68) for defin- 
ing the blue edge of the selection region. As suggested by the 
TRILEGAL model (Sect. 3.2.3 i, there probably is a bias agains 
metal-poor stars in the selection, which affects the size and shape 
of the metal-poor peak. This metallicity selection bias will be 
smoothed by the distance spread of the bulge stars. However, 
it is hard to quantify if this sharp cut-off is a real feature or 
an effect of the selection bias because of the lack of stars with 
[M/H] < -1.0 in the input MDF of the models. 



4.1 .2. Bright and faint sample stars 

It becomes clear already from Fig.|6]that there are possibly two 
peaks in the MDF: one at 0.55 and one at ~ 0.2. This double- 
peak structure becomes even more apparent if bright and faint 
stars are considered separately (Fig. |7]). The cut between bright 
and faint stars was introduced at Jq - 12™7, which should safely 
exclude RC stars from the bright sub-sample. The motivation to 
introduce a brightness cut is two-fold: first, the uncertainty in the 
metallicity is larger (~ 0.30 dex) for the coolest and brightest 
sample stars. Thus, with a brightness cut we exclude stars with 
the most uncertain metallicity determination. Second, a bright- 
ness cut is equivalent to a cut in evolutionary state, at least when 
the distance spread within the bulge of ~ 0'."7 is disregarded. 
Hence, with our sample we can check whether the MDF varies 
with evolutionary state or not. 

The most striking feature in Fig. |7] is that among the faint 
sample stars (those with Jq > 12™7) the dichotomy in metal- 
licities becomes even clearer than in the complete sample. To 
make this dichotomy more quantitative we performed a decom- 
position by Gaussian mixture models. The parameter space of 
two Gaussian mixture models was explored using a Markov 
Chain Monte Carlo simulation and a maximum-likelihood ap- 
proach to define the best-fit parameters (for details on the code 
see Nataf et al. (2011), for an in-depth explanation of MCMCs 
see ~ ' ~ ' r^^^~ _ , 

68% 



Johnson J. A. et aL| ( |201 l| l). The chosen confidence level is 



Stars with [M/H] < - 1 .4 were omitted from the fit because 
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we think these stars might be halo contaminants. The results of 
this Gaussian decomposition are presented in Fig. [8] We find 
mean metallicities of the populations of [M/H] = -0.57 + 0.03 
and [M/H] = +0.30 ± 0.04, and dispersions of 0.27 + 0.02 and 
0.28 + 0.03 dex for the metal-poor and metal-rich population, 
respectively. The metallicity difference between the two popula- 
tions, the systematically and statistically more robust figure, is 
A[M/H] = 0.87 ± 0.03. The number of stars belonging to those 
peaks is 198 + 12 and 88 + 12, respectively, yielding a fraction 
of the metal-rich population of (30.8 ± 4.2)%. The TRILEGAL 
model suggests that there could be a selection bias against metal- 



poor stars (Sect. 3.2.3 i; if true, then the fraction of metal-poor 
stars at this Galactic latitude has to be regarded as a lower hmit, 
and that of the metal-rich ones as an upper limit. 

4.1.3. Comparison with tine literature 



Recently, MDFs of bulge stars were presented by Hill et al. 
( [201 1 ) and B ensby et al.|(|2011| l, both of which find a double- 
peaked MDF. ]Hill et al. '{ 2011| l presented [Fe/H] measurements 
for 219 bulge red clump stars in Baade's window (b - -4°), 
derived also from FLAMES/GIRAFFE spectra. Their Gaussian 
decomposition yields a broader peak at [Fe/H] » -0.30 and a 
relatively narrow peak at [Fe/H] ^ -hO.32, respectively. On the 
other hand, iBensbyetaL ( 201 l| l measured from high-resolution 



UVES/VLT spectra the metallicities of micro-lensed dwarf and 
sub-giant stars scattered throughout the Galactic bulge at lati- 
tudes between b = -1° and -5°, at various longitudes. The 
dwarfs show also a clearly bi-modal distribution with peaks at 
[Fe/H] X! -0.60 and +0.3. Our result is in excellent agree- 
ment with that of Bensby et al. (2011 1, even though our field 



is several degrees (corresponding to ~ 1 kpc) away from the 
dwarfs analysed by these authors. Only the metal-poor popu- 
lation identified by Hill et al. (201 1 1 has a higher metallicity by 
~ 0.3 dex. Unfortunately, the number of analysed bulge dwarfs 
is still too low to reliably estimate the widths of the two peaks 
in their MDF In contrast to the results of Hill et al.| ( |201 1"] ), the 
peaks in the MDF found in our sample have very similar widths. 
In general, there are interesting analogies between our results 
and those of [Bensby e"taL] f20ir) and Hill et al.[ ( [20n) , and we 
conclude that there are two populations present in the GB, sepa- 
rated by 0.6 - 0.9 dex in metallicity. The small differences might 
be a result of the different selection procedures, the difference 
in evolutionary state of the samples, and the different analysis 
methods applied in the studies. 



4.2. Evolutionary effects on the MDF 

Another interesting feature of Fig. |7] is the probable lack of 
metal-rich stars among the bright sample stars (Jq < 12"7). In 
the metallicity range +0.1 < [M/H] < +0.7, eleven stars are 
bright and 72 stars are faint (13.3 + 4.3% bright stars, the un- 
certainty being calculated from Poissonian noise), whereas in 
the range -0.5 < [M/H] < +0.1 the fraction is 29.5 + 5.4% 
(39 bright and 93 faint stars). The fraction of bright stars in the 
whole sample is 22.9 + 2.7% (92 bright and 309 faint stars). This 
reduced fraction of bright stars among the metal-rich stars could 
be caused by a selection effect. We therefore inspected the frac- 
tion of bright stars in the selection region of the Galaxy models 
as a function of metallicity to check for differences in the metal- 
licity bias between faint and bright RGB stars (i.e. less and more 
evolved stars). However, no such trend is found, which suggests 
that the metallicity biases equally apply to faint and bright stars. 




Fig. 8. Metallicity distribution of the faint sample stars (Jo > 
12™7) together with the two best-fitting Gaussian mixture mod- 
els. The metallicity distribution in this brightness range is clearly 
bimodal. 



Apparently, the shape of the MDF depends on the evolution- 
ary state of the sample stars one is considering. The simplest 
explanation for this finding is that the metal-rich stars do not 
evolve all the way from the main sequence up to the RGB tip 
or the early AGB, but get "lost" somewhere on the way (e.g. 
due to enhanced mass loss; [Castellani & Castellani| [r993[l . A 
similar scenario was put forward by Chiappini et al. ( 2009[ l to 
explain the shift of 0.3 dex to lower values in the planetary neb- 
ulae oxygen abundance distribution with respect to the distribu- 
tion of giant stars. Because of the still relatively small number 
of stars in the present sample this conclusion has to be taken 
with some caution. Also, some (presumably metal-rich) stars on 
the red side of the selection region have not been observed in 
the present programme (see Fig. [TJ. However, in Sect. 4.7.3 



we 

will present another piece of evidence from the radial velocity 
dispersion of evolved stars that strongly supports the scenario of 
metallicity-dependent mass loss. 

4.3. The metallicity of the red clump stars: Is the double RC 
only connected to the metal-rich population? 

Figure [9[ shows the MDF of the bright and faint RC, with the 
definition of the two RCs as introduced in Sect. 12.21 It is clear 
from that figure that the faint RC has a much more prominent 
low-metallicity peak than the bright RC, which in turn shows 
both peaks. Their mean metallicities are [M/H] = -0.18 and 
-0.41, respectively. This would suggest that the two RCs have 
intrinsically different metallicities, in contradiction to findings 
by jde Propris et al.[ ( [201 1| , and would therefore not represent 
the same underlying population of stars. However, we believe 
that this seemingly different mean metallicity is actually caused 
by our selection procedure. Because of the way how we se- 
lected the stars along two isochrones, which become redder 
at brighter magnitudes, there are preferentially redder stars se- 
lected from the bright RC, and bluer stars from the faint RC 
(see also Fig. [5[). At a given evolutionary state, redder stars are 
more metal-rich than bluer ones. Hence, we selected preferen- 
tially more metal-rich stars in the bright RC than in the faint RC. 
This interpretation is also favoured by the RV distributions of the 
two RCs, which are practically indistinguishable (cf. Sect. 4.6 1. 
Introducing an additional colour criterion does not strengthen 
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Fig. 9. Metallicity distribution of stars belonging to the faint (up- 
per panel) and bright (lower panel) RCs, respectively. 
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Fig. 10. Luminosity functions of the metal-rich ([M/H] > -0.2, 
upper panel) and the metal-poor ([M/H] < -0.2, lower panel) 
populations around the red clump brightness. 



this conclusion because of the small number of stars that would 
be compared then. 

In fact, we suggest that the double RC feature is connected 
only to the metal-rich population. A direct hint may come from 
the luminosity functions of the metal-poor and metal-rich pop- 
ulations around the RC brightness, see Fig. 



Together with the presence of a metallicity gradient in the 



10 The metallic- 



ity cut was introduced at [M/H] = -0.2 for this diagram. The 
metal-rich population has two peaks in the Jq histogram, one at 
Jo ~ 13™2 and one at Jo ~ 13"9. It stands to reason that these 
two peaks are the manifestation of the front and back arms of the 
X-shaped bulge. The fainter peak contains fewer stars because 
its stars are at larger distance from the Galactic plane where the 
density of stars is lower ( jSaito et al.||20Ilj ). The metal-poor pop- 
ulation has one very broad, symmetric peak at Jq ~ 14™0 with 
a decline towards brighter magnitudes. Because the RC and the 
RGB bump have the same magnitude in Jo at a metallicity of 
[M/H] « -0.7 (Nataf et al., 2011), this peak may also contain 
some RGB bump stars of the metal-poor population. A statisti- 
cal t-test shows that the probability that the two luminosity func- 
tions are drawn from the same parent distribution is only 5.1%. 



Also Ness et al. (2012), from an analysis of spectra of clump 



stars along the minor bulge axis, find that a split in the luminos- 
ity function appears only for stars with [Fe/H] > -0.5, in good 
agreement with our finding. 

Another argument for the suggestion that the double RC fea- 
ture is connected only to the metal-rich population may come 
from the analogy of how the double RC feature diminishes with 
increasing distance from the Galactic plane and the presence of 
a metallicity gradient in the bulge. CMDs based on 2MASS pho- 
tometry in three different fields along the minor axis of the bulge 
(at b - -7.0°, -8.5°, and -10.0°, from top to bottom) are shown 
in Fig. 



11 



Fields with radii of 20', 30', and 40', respectively, 
were used to keep the total number of stars approximately con- 
stant. On the right hand side of each panel, a magnitude his- 
togram of the stars inside the box marked by thick lines is shown. 
The luminosity function was fitted with Gaussian and exponen- 
tial functions to reproduce the RCs, RGB bump, and the underly- 
ing distribution of RGB stars. It is clear from that figure that not 
only the peak brightness positions of the RCs are moving apart 
when going away from the plane, but also the relative number 
of stars making up the double RC feature is decreasing with dis- 
tance from the plane. 



bulge (Zoccali et al. 2008), this suggests that the double RC 
feature is connected only to the metal-rich population, which 
is much more concentrated towards the Galactic plane. In both 
Hill et aL (2011 1 and Bensby et al. (2011 1 the stars evenly dis- 
tribute among the two sub-populations. In our sample, the metal- 
rich population makes up only ~ 30% of the stars (upper limit), 
which would suggest that the fraction of the populations varies 
with Galactic latitude. Hence, the metal-rich population seems 
to be more concentrated towards the plane than the metal-poor 
one. This would also explain the metallicity gradient found by 
|Zoccali et al.| (|2008 ): If the fraction of metal-rich stars decreases 
with increasing distance from the plane, also the mean metallic- 
ity would shift to lower values. Hence, rather than a decreasing 
mean metallicity within one monolithic bulge population, proba- 
bly it is the varying percentage of two sub-populations that mim- 
ics a metallicity gradient in the bulge. 

4.4. Comparison with simulated MDFs 



Our observed MDF, corrected for sampling effects (Sect. 4.1.1 
is compared to the simulated MDFs from the Besangon and 
TRILEGAL models in Fig. [12] The simulated MDFs do not 
reproduce the observed MDF in all aspects. The MDF of the 
TRILEGAL model has clearly a much higher mean metal- 
licity than the observed mean ([M/H] - -1-0.09, compared 
to -0.34). The number of metal-poor stars is underestimated, 
whereas the super-solar metallicity stars are overrepresented in 
the TRILEGAL mo del. Clearly, the shift of -hO.3 dex app l ied by 
Van hollebeke"eraL] ( p009] ) to the MDF of Zoccali et al.| ( j2003] ) 
does not reproduce the MDF of our sample (which is at -10° 
from the Galactic plane). This probably means that an MDF cali- 
brated on fields at low Galactic latitude cannot be simply applied 
to higher latitudes because of the metallicity gradient. The mean 
metallicity in the Besangon model, on the other hand, is close to 
the observed one ([M/H] = -0.23). For this model, the distribu- 
tions of the bar (red dashed line) and the thick bulge (red dotted 
line) are shown separately in the lower panel of Fig. [12] The 
bar stars reproduce the metal-rich peak of the observed MDF 
relatively well. However, the two components are not separated 
in metallicity as much as they are in the observations. If the 
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Fig. 11. Colour-magnitude diagrams of fields at {l,b) - 
(0°,-7.0°) {top panel), (l,b) = (0°,-8.5°) (middle panel), and 
(/, b) - (0°, -10.0°) {bottom panel). For each field, a magnitude 
histogram of the stars in the thick box in the CMD is shown in 
the right hand panel, including a fit to the double RC, the RGB 
bump, and the underlying distribution of RGB stars. The thick 
horizontal bars in the magnitude histograms mark the peak mag- 
nitudes of the two RCs. 




Fig. 12. Comparison of the observed MDF corrected for sam- 
pling effects {upper panel), with the ones predicted by the 
Besanfon and TRILEGAL models in the selection region {lower 
panel). The histograms are normalised to the same area. In the 
lower panel, the red dotted line shows the MDF of "thick bulge" 
stars in the Besan9on model, whereas the dashed red line is the 
one for the bar component. 



bulge component in the Besan9on model was shifted to lower 
metallicities by 0.2 to 0.3 dex, it would very nicely reproduce 
the metal-poor peak in the observed MDF (cf. Sect. 4.7. We 
remind that these stellar population synthesis models assume 
certain metallicity distributions for each Galactic component. 
Self-consistent chemical evolution models such as lBallero et al.l 
( |2007| l are needed to do a detailed comparison. 



4.5. a-element abundances 

An important piece of evidence for the understanding of the 
GB population is the abundance of the a-elements because they 
convey information of the star formation history. In Fig. 13 we 



present the abundance of [a/Fe] as a function of the iron abun- 
dance [Fe/H]. The metal-rich and metal-poor stars also differ 
in this diagram: The metal-rich stars have clearly a lower over- 
abundance in the a-elements than the metal-poor ones. Among 
the metal-rich stars, most stars fall in the range [a/Fe] - +0.Q5 
to H-0.25, while among the metal-poor stars many stars can be 
found at over-abundances of -1-0.30 dex and more. 

Because of proximity to our field, we over-plot the abun- 
dances measured by Gonzalez et al. (2011a) in the field at 
b - -12°. To convert the individual abundances to a combined 
[a/Fe] value, we computed the weighted mean of their Ca, Ti, 
Si, and Mg abundances, using as weight the contribution to the 
absorption in the synthetic spectrum that we applied in the deter- 



mination of the ff-abundances (Sect. 3.1.3 i. The match between 



our data and those of [Gonzalez et al. ( 2011a|) is not perfect. 
The over-abundances measured by Gonzalez et al. (201 la I are 
on average slightly higher than in our data, even though a shift 
was applied to our data to take into account a systematic offset 
(Sect. 3.1.41. It is also clear that there are many more metal- 
poor stars (-1.5 < [Fe/H] < -0.5) in our sample than in that of 
[Gonzalez et al.| ( |201 la| l, and also at the metal-rich end we find 
somewhat more stars than they do. 

We also confirm the trend of decreasing a-element over- 
abundance with increasing metallicity, as was found already by 
a number of previous studies (e.g. |Fulbright et aL]|2007[|Ryde et| 
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Fig. 13. Abundance of [a/Fe] as a function of [Fe/H]. Filled 
black circles: data from this study; open red circles: data of the 
field at b - -12° from |Gonzalez et ar] ( |201 la) . The dashed line 
is a linear fit trough the data of this paper. 



aL| [20T0l|Alves-Brito et aTl[20T0l[IMerdl [2011] [Gonzalez et| 



aLl |2011a[ " ]johnson C. etaLl|2011[TBensby et al.||26lOa| i. A Hr> 
ear regression including all sample stars yields a slope of -0.12 
(dashed line in Fig. 13 i. The Icr scatter around this slope is only 
~ 0.08 dex, which suggests that it originates only from measure- 
ment error, not from real abundance scatter within the sample, at 
a given metallicity. At the high metallicity end, the trend of de- 
creasing a-element over-abundance might be shallower or even 
reversed. This seems to disagree with the data of [Gonzalez et| 
al. ( 201 la[ l overplotted in Fig. 13 however On the other hand, 
a flat trend with [Fe/H] is also found in metal-rich, microlensed 
bulge dwarf stars for some a-elements, e.g. Ca, Ti, Si, and Mg 
(Fig. 10 of Bensby et al. 201 1 1. Only oxygen is clearly found to 



continue its decrease at high metallicity. If this is a real feature 
or the result of uncertainties needs to be investigated with more 
precise line-by-line studies and larger samples, respectively. 

4.6. Radial velocities 

The RV distribution is an invaluable constraint on the nature 
of the GB, i.e. to decide whether it is a classical or a pseudo- 
bulge. The mean RV of our whole sample is -8.3 + 3.8 km s ', 
with a velocity dispersion ctrv - 76.1 ± 2.7 km s"' (Note: In 
the following, the quoted uncertainties in the mean velocity and 
the standard deviation of the distribution are always the purely 
statistical uncertainties.) The distribution of RVs of our sample 
stars is displayed in Fig. 



14 



No cold streams are obvious in 
our sample. Excluding the 50 foreground candidates identified 
in Sect. |2.3| by their high proper motion (and also the 10 stars 
for which no proper motion measurement was found), we get 
{RV} = -8.9 + 4.2kms-i and ctrv = 77.1 + 3.0kms"\ hence 
no significant change. 



We can compare our observed cr^y to those of Babusiaux et 
[ar]( [20T0) and the BRAVA team (Howard et al.l|2q08||2009)l. The 
comparison is done in Fig. [15] which is a reproduction of Fig. 7 
of Babusiaux et al. (2010 1. We emphasise that the sample se- 
lection for the BRAVA survey, Babusiaux et al. ( 2010[ l, and for 



the present study are quite different. From the BRAVA survey, 
only the results for fields along the bulge minor axis (/ a; 0°) are 
shown. Our measurement of ctrv is the lowest one in this fig- 
ure, but within the error bar of the data point at b = -12° from 
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Fig. 14. Radial velocity distribution of our sample 
graph is a Gaussian fit to the data with mean and 
ation as indicated in the figure. 



stars. The grey 
standard devi- 



140 



120 



This work 

Bobusioux et ol. 2010 
Howord et al. 2008-1-2009 
Zhao 1996 



> 100 
tr 

b 



80 



60 




Fig. 15. Comparison of observed velocity dispersions as a func- 
tion of Galactic latitude b along the bulge minor axis (/ - 0°) 
with the model of Zhao ( 1996| l. Cf. Fig. 7 of Babusiaux et al. 
( |20T0l l. 



[Babusiaux et al. (2010). Figure 15 also includes the prediction 
from the model bjTlZhao ( 1996), which is a "3D steady-state stel- 
lar dynamic model for the Galactic bar", hence a pseudo-bulge 
model. The observed ctrv as a function of b is well described 
by the model of |Zhao| ( [T996| , except for the two BRAVA data 
points farthest from the Galactic plane. The ctry of our sam- 
ple is only slightly higher than what is predicted by this model. 
Furthermore, our measurement is also in excellent agreement 
with the predictions of the pseudo-bulge model developed by 



Shen et al. (2010 their Fig. 2, lower right panel). 



The kinematics of the RC stars is an important test for the 
understanding of the nature of the double RC feature. Figure 16 
display s the RV distribution of the stars in the two RCs as defined 
in Sect.|Z2] The mean RVs are -10.4+7.3 and -6.0+6.1 km s 
and the standard deviations are 71.6 + 5.2 and 74.1 + 4.3 km s~\ 
respectively. Although the velocity dispersion of the faint RC is 
slightly larger than that of the bright RC, the two distributions 
are indistinguishable, within the error bars. To make this result 
more significant, we derived probabilities that a particular star 



13 



Uttenthaler et al.: The structure and formation of the Galactic bulge 



12 
10 
8 
6 
4 
2 

■E 

<u 

£ 10 
c 

- 8 

B 6 

o 4 

^ 



foint RC 



ru 



A 



bright RC 



Jl 



n nn 



-200 



-100 



100 



200 



RV (km/s) 



Fig. 16. Radial velocity distribution of the stars in the bright and 
faint RC, respectively. 

Table 2. Radial velocity distributions of different sub-samples. 
The top group are the values measured from the observed sam- 
ple, the bottom group are the values derived from ten simulation 
runs of the Besangon model (Robin et al. 2012 1. 



Selection 


N 


(RV) (kms-') 


ctrv (km s ' ) 


Observations 












whole sample 


401 


-8.3 H 


t3.8 


76.1 H 


t2.7 


fi < 20 mas 


341 


-8.9 H 


t4.2 


77.1 H 


t3.0 


bright RC 


96 


-10.4 H 


t7.3 


71.6H 


t5.2 


bright RC, weighted 


96 


-10.3 H 


t7.3 


71.3 H 


t5.2 


faint RC 


150 


-6.0 H 


t6.1 


74.1 H 


t4.3 


faint RC, weighted 


150 


-5.1 H 


t6.1 


73.5 H 


t4.1 


metal-rich third 


128 


-7.1 H 


t4.6 


52.3 H 


t3.3 


metal-poor third 


128 


-5.4 H 


t 8.0 


90.6 H 


t5.7 


metal-rich pop.' 


88 


-6.1 H 


t5.3 


47.7 H 


t5.6 


metal-poor pop.' 


198 


-9.1 H 


t 6.4 


90.3 H 


t6.5 


Besan^on model 












whole sample 


6933 


-11.1 H 


tO.9 


77.8 H 


tO.7 


fi < 20 mas 


6555 


-U.Oh 


t 1.0 


79.6 H 


tO.7 


bright RC 


1408 


-10.0 H 


t2.0 


74.1 H 


t 1.4 


faint RC 


2802 


-11.1 H 


t 1.5 


80.9 H 


t 1.1 


metal-rich third 


2311 


-13.1 H 


t 1.4 


65.9 H 


t 1.0 


metal-poor third 


2312 


-9.2 H 


t 1.8 


87.2 H 


t 1.3 



Notes. (1): Metal-rich and metal-poor population according to the 
Gaussian decomposition in Sect. 14.1.2) 



belongs to either the front or back over-density (bright or faint 
RC). For this end we fitted two Gaussians (for the two RCs) and 
an exponential function (for the underlying RGB stars) to the 
magnitude histogram to derive these probabilities. This proce- 
dure was applied to a 40' diameter field to increase the statistics 
for a Gaussian fitting of the two RCs. Using these membership 
probabilities as weights, we get {RV} = -10.3 + 7.3 km s"' and 
{RV} = -5.1 + 6.0kms-i, and ctrv = 71.3 ± 5.2kms-' and 
Cry = 73.5 +4. 1 km s"' for the bright and faint RC, respectively. 
Even in the weighted mean, the kinematics of the two RCs re- 
main indistinguishable. Furthermore, a statistical t-test yields a 
probability of 64.2% that the two distributions are drawn from 
the same parent population. 

The above discussed mean RVs and velocity dispersions of 
the different sub-samples are summarised in Table |2] 
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Fig. 17. Radial velocity versus metallicity. Top panel: observed 
sample stars. Foreground star candidates are represented by red 
squares, all other stars by grey squares. Bottom panel: simulated 
stars in the Besangon model (Robin et al. 2012). Blue squares 
represent the thick bulge, black squares the bar, filled red squares 
the thin disc, open red squares the thick disc, and green symbols 
represent the halo population. The over-plotted eiTor-bar sym- 
bols in both panels show the mean and standard deviation of the 
RV over the metallicity range between the vertical dotted lines. 



4.7. Combining kinematics witii abundances 

The upper panel of Fig.fTTjshows how the sample stars distribute 
in the RV vs. [M/H] plane. It it immediately clear from this dia- 
gram that the RV distribution becomes broader with decreasing 
metallicity and that at low metallicity more extreme RV values 
are found than at high metallicity. This is also illustrated by the 
large over-plotted "error-bar" symbols, which show the mean RV 
and velocity dispersion over the metallicity range between the 
vertical dotted lines. The mean RV stays more or less constant, 
whereas crgy increases monotonically with decreasing metallic- 
ity. This is further evidence that there are two different popula- 
tions present in our sample, each with its own velocity distribu- 
tion. Whereas the increase in the velocity dispersion agrees with 
the results of de Propri s et al.| ( |201 l| l, the lack of a trend of the 
mean RV with metallicity is in clear contrast to their results. This 
also holds if only the RC stars in our sample are considered. This 
is somewhat surprising because the field studied by de Propris et| 
al. (201 11) is only two degrees closer to the Galactic plane than 
our field, along the bulge minor axis. At this point we do not 
have an explanation for this discrepancy. 
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The metal-rich stars on the near side of the X (stars in the 
bright peak in the upper panel of Fig. 10 1 may be expected to 



differ in their mean radial velocity from those on the far side 



(faint peak in Fig. lOl. We find (RV) = +1.5 + 7.5 kms^' and 



ctrv = 53.4 ± 5.3 km s for the 51 bright, metal-rich RC stars 
and (RV) = -10.7 ± 8.2 km s"' and ctrv = 51.5 + 5.9 km s"' for 
the 39 faint, metal-rich RC stars, respectively (using a metallicity 
cut at [M/H] > -0.2 and the brightness criteria for the RCs of 
Sect. |2.2| i. The error bars on the mean RVs overlap, hence the 
difference is not significant. Note, however, that RGB/early AGB 
stars brighter than the RC located on the far side and stars fainter 
than the RC located on the near side will dilute any difference 
between these two groups. A still larger sample would be needed 
to confirm this difference in mean RV. 



4.7.1. Comparison with tine Besangon model 

It is interesting to compare this result with the kinematics as pro- 
vided by the simulations to see whether the BGM, based on a de- 
composition in several populations, mainly disc, bulge and bar in 
the central region, is in reasonable agreement with our data. The 
lower panel of Fig. 



17 shows for comparison the distribution of 



simulated stars in the BGM in the RV vs. metallicity plane. The 
width of the RV distribution as well as the extreme RV values 
agree well with those in the observed sample. However, the trend 
of increasing velocity dispersion with decreasing metallicity is 
much smaller in the simulation than in the observations. The 
main reason for this is that the bulge population (blue symbols 
in the lower panel of Fig. 17 1, whose velocity dispersion is rel- 
atively large, is close in metallicity to the bar population (black 
symbols), which has a smaller velocity dispersion. If the bulge 
population in the Besangon model was shifted to lower metallic- 
ities by 0.2 to 0.3 dex, it would nicely reproduce this trend, and 
also bring the range of metallicities in the model closer to the 
observed range. This is in line with the conclusions drawn from 
the comparison of the MDFs (Sect. 4.4 1. 

To compare the observed with the simulated velocity dis- 
persions, and to decrease the Poisson noise in the model, we 
calculated ten simulation runs of the Besangon model for our 
field. From these ten simulations, we obtain 6933 stars in the se- 
lection area. The results of this exercise are summarised in the 
lower part of Table |2] The mean RV as well as the velocity dis- 
persion of the simulation agree well with the observed values 
for the whole sample. Excluding stars with a total proper mo- 
tion fi > 20 mas/yr, the velocity dispersion slightly increases in 
both observations and simulations. As discussed in Sect. |3.2.3| 
this selection by proper motion is very efficient in removing the 
foreground disc. This increase in velocity dispersion is not sig- 
nificant in the observed sample, but it is significant in the sim- 
ulation. In Fig. 18 we show the histogram of the distribution of 
radial velocity from the simulations, with (dashed line) and with- 
out (solid line) the proper motion selection. As can be seen from 
this figure, the proper motion selection removes a cold (disc) 
component. 

Table |2] also lists the mean RV and dispersions of the simu- 
lated bright and faint red clump sub-samples, as well as for the 
metal-rich and metal-poor 33.3% quantile. In both the observa- 
tions and in the simulations the bright RC has a slightly smaller 
velocity dispersion; while this increase is again not significant 
in the observed sample, it is significant in the simulations. The 
reason for this probably is that, as already discussed in Sect. |4.3| 
the faint RC is more contaminated by stars from the metal-poor 
population, which has a higher velocity dispersion than the bar 
population, the probable carrier of the double clump feature. 
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Fig. 18. Radial velocity distribution of stars in ten simulation 
runs of the Besangon model ( [Robin et al.j |2012| l for our field. 
The full line is for all stars in the selection region, the dashed 
line shows the distribution for stars with total proper motion 
fi <20 mas/yr 



The metal-rich sample has a smaller velocity dispersion than 
the metal-poor one. This is due to the fact that the metal-rich 
one is dominated by the bar population, while the more metal- 
poor one is dominated by the bulge. The velocity dispersion of 
the metal-poor third of the simulation agrees nicely with the ob- 
served one (87.2 km s"' to be compared with 90.6 km s"', see 
below). However, the metal-rich third has a larger dispersion 
in the simulations (65.9 km s"') than in the observed sample 

the metal- 
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(52.3 km s"'). As seen in the lower panel of Fig. 
rich part of the sample is dominated by the bar, which follows 
Fux ( 1999) dynamics. But there is also a contribution from the 



disc, which is mainly removed using the proper motion selec- 
tion. In the very low metallicity tail, there are a few halo stars 
and a noticeable contribution from the thick disc. It has, in our 
model, a smaller dispersion than the bulge, which has about the 
same metallicity. The velocity dispersion of the thick disc in the 
inner galaxy is not known and the simulation we are proposing 
here is just an attempt to extrapolate the thick disc as seen in the 
solar neighbourhood to the central region. Though, the disper- 
sion we obtain for this intermediate and low metallicity bin is in 
very good agreement with the data. 

We conclude that this new model, having both a flaring bar 
and a bulge component in the central region, explains very well 
the observed relation between the metallicity and the kinematics 
in this bulge field. 

4.7.2. Comparison witli tine Fux model 

The increase of ctrv with decreasing metallicity is reminiscent 
of the results of |Babusiaux et al!]p010 ), who also found a clear 
distinction in kinematics between the metal-poor and the metal- 
rich stars in bulge fields at different Galactic latitudes. In a re- 
production of their Fig. 8, we show in Fig. 19 a comparison 



with the bulge dynamics as predicted by the 3D self-consistent 
N-body baiTed models of the Milky Way by |Fux| ( |1999| l. As 



in [Babusiaux et al. (2010 



ctrv is plotted in different sym- 
bols for the complete sample (diamonds), the 33.3% metal-poor 
quantile (open downward-pointing triangles), and the 33.3% 
metal-rich quantile (filled upward-pointing triangles). We find 
ctrv - 52.3 ± 3.3 km s"' for the metal-rich third of the sample 
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Fig. 19. Observed RV dispersions as a function of Galactic lat- 
itude b along the bulge minor axis, compared to the models 
of |Fux] ( [T999l , cf. Fig. 8 of 'Babusiaux et al.|(pOTO) The blue 
symbols are the data points from Babusiaux et al. ( 2010| l, the 
black symbols are from this work. The diamond symbol repre- 
sents the velocity dispersion of the complete sample, while the 
open downward-pointing triangle represents the 33.3% metal- 
poor quantile and the filled upward-pointing triangle represents 
the 33.3% metal-rich quantile, respectively. 



and (Trv - 90.6 + 5.7 km s ' for the metal -poor third (Table H. 
By adding only one more data point at ^ = -10° we can put 
much more weight on the conclusions of Babusiaux et al. ( 2010| l. 
Within the eiTor bars, ctrv of the metal-poor quantile is constant 
throughout the bulge, from b = -4° to -12°. The spheroid model 
of pux ( 1999 ) also predicts crgy to be fairly constant as a func- 
tion of b, except for a dip around b = -9°. However, the ob- 
served (Trv is systematically lower than the ctrv predicted by 
the spheroid model. The velocity dispersion of the 33.3% metal- 
rich quantile, on the other hand, decreases monotonically with 
increasing distance from the Galactic plane. This follows closely 
the prediction by the disc/bar model of Fux ( 1999 ), but except at 
b - -4° it is also always below the predicted value. 

We must caution at this point that, while the definition of 
quantiles is a robust measure, the metallicity cuts will be at 
different absolute values, depending on the metallicity scale 
used in a given study and on the metallicity distribution in 
the different fields. We therefore also derived the RV disper- 
sions of the two populations from the Gaussian decomposi- 
tion done in Sect. 



4.1.2 



Assuming that both have Gaussian 
RV distributions, we find that the metal-poor population has 
(RV) = -9.1 + 6.4kms-i with ctrv = 90.3 + 6.5 kms"', while 
the metal-rich population has (RV) = -6.1 ± 5.3 km s ' and 
CTRV = 47.7 + 5.6kms-i (Table|2]i. 

4.7.3. Evidence for metallicity-dependent mass loss from the 
velocity dispersion of evolved stars 

In Sect. |4.2| we speculated that the lack of metal-rich stars among 
the bright sample stars might be caused by strong, metallicity- 
dependent mass loss that preferentially removes metal-rich stars 
from the evolutionary pathes (the so-called AGB manque). There 
is another, even stronger piece of evidence coming from the 
RV distributions of AGB stars and planetary nebulae (PNe) 
in the PG3 field that supports this explanation. RVs of Mira 
and semi-regular variables (SRVs) in the PG3 have been pre- 
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Fig. 20. Velocity dispersion of the metal-rich and metal-poor 
populations (Table |2|, as well as of several samples of evolved 
stars towards the PG3 field: AGB stars from Uttenthale r et al.l 
(2007 U07), Miras and SRVs from Schultheis (1998) without 
and with selection using a period - /T-magnitude relation, and 
PNe from Durand et al.| ( |199 8 ). The grey shaded areas mark the 
range of uncertainty of the metal-rich and metal-poor popula- 
tion, respectively. 



sented by |Schultheis| ( [T998|) and [Uttenthaler etaL] p)07] l. The 
27 AGB stars studied by [Uttenthaler et al.| ( |2007| l by means of 
high-resolution UVES /VLT spectra have a velocity dispersion 
of 90.9 + 12.6 km s"'. [Schultheis (1998 ) found ctrv = 95.2 ± 
8.8 km s-i for 59 PG3 Miras and ctrv = 83.0 ± 6. 1 km s"' for 94 
SRVs. If additional criteria based on the period - /iT-magnitude 
relation are applied to select only probable bulge members, these 
numbers change to ctrv = 110.4 ± 13.2 km s"' (Miras, = 36) 
and o-Rv = 89.4 + 9.6 km s^' (SRVs, = 43). Furthermore, 
we retrieved RV measurements of PNe from the compilation of 
Durand et al.[ ( |1998| . Nineteen PNe were found within a search 
radius of 4° from the centre of our FLAMES field. Those 19 PNe 
have a velocity dispersion of 94.1 + 15.7 km s"'. Although this 
last number is not highly significant owing to the small number 
of objects, it is clear that these velocity dispersions agree very 
well with the one determined for the metal-poor sub-population, 
but disagree with that of the metal-rich one. This statement is 
summarised in Fig. |20] We conclude from this that the AGB 
stars and PNe in the PG3 field descend mainly from the metal- 
poor sub-population. Or, to say it the other way round, we con- 
clude that the metal-rich stars probably never evolve up to the 
AGB stage and beyond because they terminate their evolution 
before. This conclusion has also important consequences on the 
interpretation of the dual chemistry found in many bulge PNe 
(Guzman-Ramirez et al. 201 1\ , a phenomenon which was pre- 
viously assumed to be a consequence of high-metallicity PNe 
precursor stars ( Perea-Calderon et al.(|2009|l. 

Apparently, there is a significant metal-rich population 
present out to b > -10°, however it seems to be unable to pro- 
duce highly evolved stars. Also the dearth of super-solar metal- 
licity stars among M giants in fields close to the Galactic plane 
^Rich et al. 2007 , 201 2 ) may be caused by enhanced mass loss 
off metal-rich stars, in line with our result and interpretation. 
Finally, enhanced mass loss might also be required to under- 



stand the binary fraction of low-mass white dwarfs (Brown et 
[aL][201 1| . A puzzle in this context is the presence of technetium- 
rich, long-period Mira stars towards the PG3 field ( Uttenthaler] 



et al. 2007 1, which suggest the presence of a younger and more 
massive (M x I.5M0) population in the GB, but at the same 
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time carbon stars are extremely rare or even absent in the bulge 
([Blanco & Terndrup', 19891 



Ng 1997 



SchultheisI [T998] i. This 



seems to be incompatible with the hypothesis of bulge AGB stars 
descending only from the old, metal-poor population. A solution 
to this problem could be that the Tc-rich AGB stars are the re- 
sult of blue straggler evolution. Also problematic in this respect 
are the various selection effects that are certainly present on the 
AGB and PN samples in the GB, and only more detailed inves- 
tigations, e.g. on the metal content of bulge AGB stars, can give 
more reliable answers. 

A natural consequence of the loss of metal-rich stars by en- 
hanced stellar winds would be that samples of bright M-type 
giants, as used by the BRAVA survey to probe the bulge velocity 
dispersion, could be biased towards metal-poor stars. In Sect. 4.6 
we found that the observed ctrv as a function of b is well de- 
scribed by the model of Zhao ( 1996 ), except for the two BRAVA 
data points farthest from the Galactic plane (Fig. 15 1. A bias 



towards low metallicity could be the reason why the ctrv mea- 
sured by BRAVA at large distances from the plane are systemat- 
ically higher than the predictions by the |Zhao ( 19961 model: at 
this Galactic latitude, the velocity dispersions of the metal-rich 
and the metal-poor populations differ considerably, thus a bias 
towards more metal-poor stars will yield an enhanced velocity 
dispersion. 

5. Comparison of observed and simulated CIVIDs 

Comparisons with observed CMDs are an important test of 
Galaxy models. The 2MASS all-sky catalogue is an ideal source 
for such constraints (see e.g. Robin et al. 2012jl. In Fig. 2T|we 



provide a comparison for our field. No extinction correction was 
applied to the observed CMD. Rather, we applied the Marshall 
[etaL] (|2006^ 3D extinction model to both the TRILEGAL and 
BGM simulations. In addition, photometric errors in each filter 
of the order of those in the 2MASS catalogue were added to both 
the Besangon and the TRILEGAL simulation. 

The CMD in this area of the sky is dominated by three se- 
quences: a sequence on the blue side formed by main-sequence 
disc stars in the foreground, a vertical sequence formed by red 
clump stars distributing at various distances from the sun, and 
a red sequence formed by bulge stars, along which our selec- 
tion of spectroscopic targets was made. These three sequences 
are most clearly resolved in the Besancon model (middle panel 
of Fig. 21 1, even though photometric errors of the same order of 
magnitude as in the 2MASS catalogue have been added to the 
simulated stars. 

The two red clumps in the observed CMD at 7 ~ 13™3 and 
J ~ 14™0 are not well reproduced in the simulated CMDs for 
this field. This is not surprising for the TRILEGAL model, which 
uses a triaxial structure to describe the bulge. On the other hand, 
it has been shown that the flaring bar model developed by |Robin 
|et al.| ( [20T2| can reproduce the double clump feature, at least for 
a field at (I, b) - (0°, -7°). The double clump feature is not obvi- 
ous in this single simulation plotted here, but an inspection of the 
ten combined simulations (Sect.[47T|, which were calculated to 
reduce the Poisson noise, the two clumps are well discernible at 
J magnitudes of ~ 13™6 and ~ 14™1. Hence, while the faint RC 
is well matched by the Besancon model, the bright one is slightly 
underestimated in its brightness (e.g. due to an overestimation of 
its distance). 

In the TRILEGAL CMD a vertical sequence of RC stars is 
dominating in the bulge, which is probably a result of the triax- 
ial bulge structure seen almost end-on, such that the bulge RC 
stars distribute over a relatively large range of distances from 



the sun. Furthermore, the number of stars in the disc sequence 
seems to be lower in the TRILEGAL simulation than it is in both 
the observed and the Besangon CMD. The number of stars in the 
CMDs of Fig. |2T|is not normalised; star counts constitute an im- 
portant test of population synthesis models ( Girardi et al. |2005[ l. 
In the colour and magnitude range of Fig. 21 the observed CMD 



contains 2148 stars, while the Besangon and TRILEGAL models 
contain 2456 and 1476 stars, respectively. 

While some details in the simulated CMDs somewhat devi- 
ate from observed CMDs, it is clear that the Galaxy models have 
reached already a very high level of accuracy and can reproduce 
most of the stellar populations at least qualitatively. 

6. Conclusions 

We presented an analysis of spectra of ~ 400 red giant stars to- 
wards a field at {I, b) - (0°, -10°). From these spectra, we derive 
metallicities, iron and a-element abundances, and radial veloci- 
ties. It is the first study that presents a homogeneous analysis of 
stars on the upper and lower RGB of the Galactic bulge. These 
data, as well as photometric data from the 2MASS survey, are 
compared to predictions by the TRILEGAL ( Girardi et al. 2005| l 
and Besangon ( Robin et al. 2012) models of the Galaxy. These 
models are also used extensively to interpret our observations. 

The mean metallicity of the whole sample is [M/H] = 
-0.34, and the radial velocity dispersion is ctrv ~ 76kms"'. 
In this study we confirm the presence of two sub-populations 
in the bulge at peak metallicities that differ by ~ 0.9 dex. The 
peak metallicities are found at [M/H] ~ -0.6 and ~ -hO.3, 
with roughly equal dispersions. The sub-populations have sig- 
nificantly different kinematics and (^-element abundances: the 
metal-rich population has a narrow velocity distribution (ctrv ~ 
50kms ') and low cc-abundances, whereas the metal -poor one 
has a broad distribution in radial velocities (ctrv ~ 90kms"') 
and high c-abundances. The metal -rich population makes up 
~ 30% of our sample. This confirms and fosters recent results 
by [Babusiaux efaLl ( |20T0| , [Bensby et al] (|20TT), and [Hill et 



al.'(| 2011[ l in most^spects. Also the kinematic models of [Zhao 
(1996)1 and Fux ( 1999) are confronted with our data, and good 
agreement is found. 

Furthermore, we find support for the suggestion by 
[ Babusiaux et al.| ( |2010| l that the ratio in number of these two sub- 
populations is a function of angular distance from the Galactic 
plane, which explains the observed metallicity gradient. The new 
scheme in the Besanfon Galaxy model explains well the fea- 
tures seen in this sample, with two populations, a flaring metal- 
rich bar and a more metal-poor classical or thick bulge, even 
though the mean metallicity of this thick bulge should be ad- 
justed downward by 0.2 to 0.3 dex. The metallicity distribution 
simulated with the TRILEGAL model yields a clearly too high 
mean metallicity to fit our observations. 

This study also presents the first medium-resolution, high- 
S/N spectra of stars belonging to the two red clumps of the bulge. 
A small difference in mean metallicity between the two RCs is 
attributed to selection effects. We find that the double RC might 
be entirely due to the metal-rich sub-population. The kinematic 
difference between the two RCs when measured with radial ve- 
locities alone is at most very small. Some difference in the mean 
radial velocity of bright and faint metal-rich RC stars is found, 
although a still larger sample would be needed to confirm that 
stars on the far side preferentially move faster towards us than 
the ones on the near side. 

We also find indications that the metallicity distribution 
function of the bulge might depend on the evolutionary state of 
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Fig. 21. Observed and simulated CMDs plotted as density maps. Left panel: observed CMD from 2MASS. Middle panel: simulated 
CMD from the Besangon model. Right panel: simulated CMD from the TRILEGAL model. The number of stars is not normalised. 
The scale on the right hand side illustrates to what colours the density of stars are connected in the diagrams. 



the considered sample stars. In particular, there are fewer metal- 
rich stars present among the brighter, more evolved stars in our 
sample. Although this result could be a consequence of selec- 
tion bias or the still limited number of stars in our sample, we 
interpret this as the result of strong mass loss that causes the 
most metal-rich stars to terminate their evolution before the most 
advanced states. A very robust indication for that comes from 
the fact that AGB stars and PN in the outer bulge have a veloc- 
ity dispersion that agrees with that of the metal-poor population 
identified here, but clearly disagrees with that of the metal-rich 
population. We conclude, hence, that metal-rich stars might lose 
so much mass that they skip phases of post-main-sequence stel- 
lar evolution. This means that the selection of M giants in the 
BRAVA survey could be biased towards metal-poor stars. We 
also confirm the trend of decreasing c-element over-abundance 
with increasing iron abundance, as found in previous studies 
(e.g. [Gonzalez et aL 2011aj . 

Our study also provides for an explanation how the existence 
of a dual bulge identified here and elsewhere can be reconciled 
with kinematic studies that conclude that the Milky Way bulge 
might be a pure pseudo-bulge (e.g. Shen et al. , 2010). In our ver- 
sion of the Besan9on model, the (thick or classical) bulge popu- 
lation has only 4% the mass of the bar, even when its contribu- 
tion is significant when one moves away from the plane. Thus it 
is compatible with the result of Shen et al. (2010|), which finds 
that the fraction of the classical bulge cannot be more than 8% 
of the disk mass. The pseudo-bulge dominates, but there is some 
space for a classical bulge, which we see in the data. 
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Table 3. Parameters and abundances of the sample stars. 



No. 


RA(J2000) 


Dec(J2000) 


J 


H 


K 


Teff 


RV 


[M/H] 


[Fe/H] 


[a/Fe] 


FG 




deg. 


deg. 


mag 


mag 


mag 


K 


kms~' 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



001 


276.814513 


-33 


.589966 


9.495 


8.524 


8.256 


3465 


+117.328 


+0.94 










* 


002 


276.544710 


-33, 


.787083 


9.512 


8.624 


8.369 


3697 


+20.277 


+0.21 












003 


276.608107 


-33 


.667065 


9.516 


8.615 


8.298 


3520 


+76.484 


-1.00 












005 


276.918684 


-33 


.751503 


9.696 


8.820 


8.527 


3629 


+55.090 


-1.12 


-0, 


.95 


+0, 


,07 




006 


276.432104 


-33 


.784958 


9.926 


9.061 


8.704 


3514 


-29.587 


+0.82 












007 


276.714893 


-33 


.891556 


9.978 


9.096 


8.838 


3698 


+ 12.913 


-0.56 


-0, 


.48 


+0, 


,11 




008 


276.600416 


-33, 


.681347 


10.007 


9.121 


8.891 


3756 


-135.108 


+0.06 


+0, 


.05 


+0, 


,32 




009 


276.807769 


-33, 


.757442 


10.057 


9.178 


8.941 


3749 


+106.494 


-1.10 


-0, 


.81 


+0, 


,01 




010 


276.873765 


-33 


.726830 


10.098 


9.183 


8.932 


3637 


-128.768 


-0.39 


-0, 


.41 


+0, 


,60 




Oil 


276.817925 


-33 


.649303 


10.115 


9.191 


8.923 


3578 


-59.147 


-0.76 












012 


276.611791 


-33 


.631435 


10.170 


9.330 


9.094 


3843 


+74.703 


-0.13 


+0, 


,06 


+0, 


,09 




013 


276.506457 


-33, 


.723846 


10.213 


9.322 


9.062 


3681 


+27.664 


-0.06 


+0, 


.02 


+0, 


,29 




015 


276.488172 


-33 


.844830 


10.232 


9.400 


9.158 


3847 


+85.083 


-0.39 


-0, 


.30 


+0, 


,05 




016 


276.682337 


-33 


.955391 


10.266 


9.421 


9.194 


3845 


+25.704 


-0.27 


-0, 


.06 


+0, 


,08 




017 


276.595604 


-33 


.694199 


10.406 


9.574 


9.312 


3803 


+50.134 


-0.54 


-0, 


.32 


+0, 


,09 




018 


276.592533 


-33, 


.625336 


10.409 


9.579 


9.366 


3914 


-59.013 


-0.70 


-1, 


.07 


+0, 


,55 




019 


276.586978 


-33, 


.867764 


10.470 


9.639 


9.438 


3934 


+178.014 


-0.85 


-1, 


.15 


+0, 


,32 




020 


276.704314 


-33 


.794868 


10.542 


9.678 


9.446 


3794 


+33.422 


+0.29 










* 


021 


276.653399 


-33 


.955818 


10.642 


9.747 


9.556 


3815 


+72.019 


-0.56 


-0, 


.31 


-0, 


,04 




022 


276.471245 


-33, 


.828197 


10.657 


9.769 


9.520 


3711 


-80.548 


+0.28 










* 


023 


276.565931 


-33 


.660713 


10.690 


9.821 


9.542 


3687 


-42.688 


+0.12 












024 


276.667873 


-33 


.869152 


10.726 


9.949 


9.659 


3857 


-7.078 


-0.14 


+0, 


.12 


+0, 


,08 




025 


276.771326 


-33 


.748791 


10.735 


9.924 


9.721 


3969 


-131.548 


-0.62 


-0, 


.29 


+0, 


,12 




026 


276.566599 


-33 


.576355 


10.747 


9.916 


9.673 


3849 


+ 15.044 


+0.02 


+0, 


.12 


+0, 


,17 




027 


276.621730 


-33, 


.925117 


10.779 


9.943 


9.708 


3849 


-0.464 


+0.31 


+0, 


.38 








028 


276.566492 


-33, 


.894100 


10.807 


10.005 


9.774 


3932 


-78.865 


-0.55 


-0, 


.19 


+0, 


,12 




029 


276.632910 


-33 


.850189 


10.843 


10.032 


9.847 


4008 


+305.574 


-0.58 


-0, 


.87 


+0, 


,23 




030 


276.876718 


-33 


.749680 


10.911 


10.121 


9.909 


3991 


+13.784 


-0.70 


-0, 


,93 


+0. 


,18 




031 


276.833702 


-33, 


.646965 


10.933 


10.095 


9.844 


3808 


-12.983 


-0.68 


-0, 


,33 


+0, 


,05 




032 


276.664170 


-33, 


.889168 


11.002 


10.132 


9.895 


3770 


+59.999 


+0.20 












033 


276.608510 


-33, 


.775723 


11.064 


10.240 


10.041 


3953 


+46.289 


-0.25 


-0, 


,02 


-0, 


,15 




034 


276.530697 


-33, 


.808578 


11.077 


10.243 


10.056 


3959 


+72.605 


-0.58 


-0, 


,60 


+0, 


,25 




035 


276.734398 


-33, 


.905167 


11.112 


10.288 


10.008 


3775 


+37.840 


-0.17 


-0, 


,06 


+0, 


,16 




036 


276.693013 


-33 


.893147 


11.119 


10.247 


10.060 


3874 


+ 110.060 


-0.54 


-0, 


,34 


+0, 


,04 




037 


276.455229 


-33 


.789265 


11.126 


10.329 


10.132 


4017 


+25.919 


-0.06 


-0, 


,44 


+0, 


,29 




038 


276.531016 


-33 


.587563 


11.157 


10.323 


10.122 


3933 


-105.342 


-0.38 


-0, 


,03 


+0, 


,22 




039 


276.476724 


-33, 


.652061 


11.191 


10.461 


10.223 


4074 


-99.668 


+0.03 


-0, 


,15 


+0, 


,13 




040 


276.632713 


-33 


.608341 


11.204 


10.397 


10.238 


4075 


-213.279 


-0.82 


-1, 


,06 


+0, 


,38 




041 


276.611649 


-33, 


.623035 


11.212 


10.433 


10.228 


4037 


+7.052 


-0.40 


-0, 


,55 


+0, 


,20 




042 


276.514543 


-33, 


.867790 


11.212 


10.396 


10.256 


4096 


+3.906 


-0.65 


-0, 


,80 


+0, 


,27 




043 


276.545896 


-33, 


.910049 


11.236 


10.429 


10.287 


4109 


-60.971 


-0.41 


-0, 


,55 


+0, 


,20 




044 


276.501519 


-33, 


.712631 


11.270 


10.456 


10.228 


3917 


-66.401 


-0.08 


+0, 


,18 


+0, 


,10 




045 


276.497877 


-33, 


.881786 


11.378 


10.535 


10.362 


3969 


+240.911 


-1.11 


-1, 


,12 


+0, 


,19 




046 


276.469519 


-33 


.651062 


11.403 


10.619 


10.362 


3921 


-1.719 


-0.54 


-0, 


,34 


+0, 


,16 




047 


276.718276 


-33, 


.890141 


11.407 


10.593 


10.388 


3957 


+29.551 


+0.05 












048 


276.593535 


-33, 


.598366 


11.415 


10.584 


10.345 


3857 


-9.884 


+0.14 


+0, 


,23 


+0, 


,33 




049 


276.751404 


-33 


.686954 


11.428 


10.661 


10.459 


4065 


+81.687 


-0.06 


-0, 


,37 


+0, 


,21 




050 


276.668148 


-33 


.844818 


11.433 


10.578 


10.353 


3829 


+29.158 


-0.44 


-0, 


,12 


-0, 


,10 




051 


276.897946 


-33, 


.802887 


11.446 


10.612 


10.445 


3992 


-53.729 


-0.83 


-1, 


,01 


+0, 


,26 




052 


276.827522 


-33 


.672203 


11.457 


10.650 


10.435 


3952 


-28.859 


+0.00 












053 


276.598353 


-33 


.763283 


11.464 


10.658 


10.412 


3893 


+ 16.198 


-0.02 


+0, 


,20 


+0, 


,09 




054 


276.753347 


-33 


.748871 


11.480 


10.655 


10.504 


4049 


+21.622 


-0.74 


-0, 


,92 


+0, 


,23 




058 


276.763338 


-33 


.726757 


11.555 


10.775 


10.634 


4167 


-133.305 


-0.97 


-1, 


,07 


+0, 


,26 




060 


276.814847 


-33, 


.799526 


11.592 


10.814 


10.653 


4126 


+ 19.204 


-0.64 


-0, 


,79 


+0, 


,18 




061 


276.745468 


-33, 


.614552 


11.610 


10.791 


10.615 


4011 


-206.421 


-0.56 


-0, 


,75 


+0, 


,24 




065 


276.767964 


-33 


.723507 


11.682 


10.895 


10.664 


3961 


-59.417 


+0.60 












067 


276.764757 


-33 


.714439 


11.695 


10.914 


10.713 


4036 


-42.029 


-0.33 


-0, 


.38 


+0, 


,13 




069 


276.732436 


-33, 


.691238 


11.750 


10.953 


10.767 


4035 


+28.948 


-0.30 


-0, 


,43 


+0, 


,20 




070 


276.809016 


-33 


.594540 


11.763 


10.955 


10.743 


3957 


-54.899 


-0.19 


-0, 


,20 


+0, 


,00 




072 


276.594486 


-33 


.730743 


11.791 


11.111 


10.937 


4325 


-52.206 


-0.52 


-0, 


,59 


+0, 


,20 




073 


276.836739 


-33, 


.717419 


11.816 


11.071 


10.883 


4139 


-56.535 


-0.28 


-0, 


,39 


+0, 


,16 




075 


276.910260 


-33 


.767544 


11.837 


11.034 


10.873 


4070 


+45.375 


-0.45 


-0, 


,59 


+0, 


,16 


* 


080 


276.663676 


-33 


.680828 


11.956 


11.185 


11.002 


4099 


+67.611 


-0.63 


-0, 


,80 


+0, 


,21 




081 


276.598738 


-33 


.673889 


11.964 


11.282 


11.096 


4294 


-155.581 


-0.12 


-0, 


,25 


+0, 


,23 




084 


276.723840 


—33 


.562702 


11.989 


11.339 


11.120 


4290 


+29.832 


+0.49 


+0, 


,21 


+0, 


,16 
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Table 3. Continued. 



No. 


RA(J2000) 


Dec(J2000) 


J 


H 


K 




RV 


[M/H] 


[Fe/H] 


[a/Fe] 


FG 




deg. 


deg. 


mag 


mag 


mag 


K 


km 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



085 
086 
088 
089 
096 
101 
102 
104 
107 
108 
110 
112 
115 
116 
117 
118 
121 
123 
128 
131 
132 
133 
136 
139 
142 
143 
144 
147 
148 
154 
156 
158 
159 
167 
169 
172 
178 
180 
182 
183 
184 
186 
187 
188 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 



276.663499 
276.810910 
276.843442 
276.708411 

276.728929 
276.650041 
276.603598 
276.764294 
276.796933 
276.818554 
276.831116 
276.844085 
276.876631 
276.886129 
276.839217 
276.605428 
276.850407 
276.733228 
276.674814 
276.862902 
276.529376 
276.754996 
276.670396 
276.669257 
276.592785 
276.720586 
276.724865 
276.605327 
276.808190 
276.884497 
276.802245 
276.607949 
276.506642 
276.696311 
276.585971 
276.674085 
276.673132 
276.715472 
276.772886 
276.791120 
276.841313 
276.721044 
276.533345 
276.753476 
276.538801 
276.627097 
276.790756 
276.717695 
276.864580 
276.694217 
276.841761 
276.749744 
276.638867 
276.601595 
276.910018 
276.655418 
276.600953 
276.718361 
276.651481 
276.778831 
276.529719 
276.854131 
276.628729 
276.798530 
276.836967 



-33.799038 
-33.617214 
-33.804676 
-33.645660 

-33.625156 
-33.653530 
-33.725780 
-33.611069 
-33.752850 
-33.721001 
-33.661152 
-33.719276 
-33.669006 
-33.817890 
-33.735954 
-33.607613 
-33.727764 
-33.686802 
-33.794342 
-33.842964 
-33.635666 
-33.640835 
-33.678337 
-33.560745 
-33.609268 
-33.791100 
-33.593685 
-33.599644 
-33.749828 
-33.711712 
-33.615543 
-33.626446 
-33.859463 
-33.687943 
-33.583626 
-33.722809 
-33.740860 
-33.633411 
-33.712326 
-33.629635 
-33.728474 
-33.780018 
-33.873089 
-33.945801 
-33.644936 
-33.841202 
-33.846645 
-33.835640 
-33.723576 
-33.601109 
-33.689163 
-33.856068 
-33.930489 
-33.643906 
-33.774948 
-33.667728 
-33.920620 
-33.917198 
-33.734238 
-33.720421 
-33.594994 
-33.731907 
-33.743752 
-33.662918 
-33.769405 



12.006 
12.009 
12.058 
12.060 
12.172 
12.257 
12.263 
12.270 
12.304 
12.315 
12.343 
12.386 
12.398 
12.401 
12.407 
12.438 
12.470 
12.482 
12.520 
12.596 
12.607 
12.630 
12.694 
12.724 
12.796 
12.801 
12.806 
12.834 
12.835 
12.903 
12.915 
12.936 
12.947 
12.992 
13.016 
13.036 
13.069 
13.078 
13.085 
13.088 
13.093 
13.110 
13.112 
13.119 
13.134 
13.146 
13.151 
13.153 
13.154 
13.163 
13.167 
13.171 
13.179 
13.190 
13.190 
13.192 
13.193 
13.194 
13.196 
13.203 
13.210 
13.221 
13.226 
13.233 
13.243 



11.199 
11.269 
11.340 
11.331 
11.457 
11.553 
11.594 
11.528 
11.596 
11.546 
11.625 
11.681 
11.613 
11.716 
11.626 
11.750 
11.850 
11.722 
11.823 
11.871 
11.953 
11.933 
11.991 
12.024 
12.171 
12.219 
12.156 
12.130 
12.127 
12.231 
12.209 
12.308 
12.320 
12.354 
12.356 
12.399 
12.451 
12.449 
12.412 
12.478 
12.437 
12.489 
12.435 
12.543 
12.484 
12.522 
12.493 
12.548 
12.547 
12.508 
12.539 
12.572 
12.505 
12.594 
12.585 
12.617 
12.561 
12.505 
12.496 
12.613 
12.613 
12.595 
12.621 
12.540 
12.643 



11.069 
11.096 
11.215 
11.179 
11.252 
11.355 
11.421 
11.326 
11.481 
11.354 
11.495 
11.584 
11.454 
11.560 
11.535 
11.607 
11.667 
11.609 
11.717 
11.725 
11.827 
11.784 
11.848 
11.849 
11.994 
11.977 
12.024 
12.014 
11.988 
12.130 
12.115 
12.179 
12.168 
12.170 
12.206 
12.232 
12.242 
12.266 
12.316 
12.337 
12.271 
12.371 
12.314 
12.370 
12.360 
12.389 
12.363 
12.410 
12.417 
12.393 
12.385 
12.406 
12.351 
12.394 
12.420 
12.458 
12.403 
12.362 
12.354 
12.452 
12.466 
12.447 
12.494 
12.415 
12.456 



4134 
4185 
4343 
4261 
4171 
4213 
4352 
4119 
4392 
4080 
4334 
4444 
4115 
4346 
4277 
4380 
4441 
4278 
4446 
4277 
4514 
4341 
4342 
4278 
4454 
4391 
4504 
4407 
4335 
4521 
4453 
4574 
4514 
4398 
4434 
4444 
4385 
4424 
4535 
4584 
4393 
4618 
4462 
4585 
4530 
4569 
4481 
4605 
4620 
4537 
4499 
4544 
4381 
4469 
4527 
4637 
4480 
4370 
4351 
4583 
4613 
4519 
4642 
4405 
4484 



+58.329 
+9.035 
+84.524 
-58.044 
+38.482 

-123.588 
-37.926 
+23.454 
+ 16.073 
-55.377 

-109.396 
-3.301 

-172.477 
-54.145 

-141.307 
+31.103 
-2.194 
-85.450 
+32.514 
-28.810 

-159.866 

+136.734 
-16.510 
-8.143 

-120.627 
-28.694 
+47.662 
-13.128 
-40.138 
-74.863 
+5.379 
+47.731 
-67.344 
+ 16.036 
-21.109 
+62.310 
-14.891 
+57.926 
-44.163 
-74.273 
+13.471 
+29.773 

-142.788 
+ 19.007 
+ 11.192 

-175.576 
-8.106 
+56.213 
-33.430 
-3.184 
-44.337 
-22.102 
+89.036 
+55.281 
-23.265 
-12.106 
-72.007 

-112.531 
-44.013 
-66.094 
+26.825 
+62.408 
-39.333 
+95.619 
+21.124 



-0.41 
+0.60 
-0.42 
-0.14 
-0.48 
-0.24 
-0.61 
+0.06 
-1.03 
+0.16 
-0.97 
-0.69 
-0.36 
-0.99 
-0.92 
-0.99 
-0.10 
-0.76 
-0.70 
-0.37 
-0.74 
-0.22 
-0.50 
+0.06 



-0.51 
+0.22 
-0.54 
-0.30 
-0.69 
-0.34 
-0.74 
-0.21 
-1.17 
-0.05 
-1.11 
-0.89 
-0.51 
-1.08 
-1.00 
-1.18 
-0.32 
-0.82 
-0.81 
-0.55 
-1.06 
-0.46 
-0.50 
-0.16 



+0.15 
+0.23 
+0.12 
+0.21 
+0.21 
+0.14 
+0.25 
+0.19 
+0.35 
+0.12 
+0.32 
+0.32 
+0.29 
+0.24 
+0.27 
+0.38 
+0.28 
+0.23 
+0.20 
+0.22 
+0.49 
+0.20 
+0.21 
+0.30 



-0.14 


-0.36 


+0.08 


-0.00 


-0.14 


+0.09 


+0.17 


-0.04 


+0.20 


-0.07 


-0.12 


+0.20 


+0.55 


+0.32 


+0.18 


-0.27 


-0.64 


+0.16 


+0.02 


-0.37 


+0.28 


+0.06 


-0.09 


+0.18 


+0.46 


+0.23 


+0.18 


+0.65 


+0.39 


+0.19 


-0.87 


-1.00 


+0.30 


+0.39 


+0.18 


+0.17 


+0.64 


+0.37 


+0.15 


-0.70 


-0.80 


+0.24 


+0.34 


+0.12 


+0.23 


-0.20 


-0.41 


+0.18 


+0.47 


+0.23 


+0.24 


-1.05 


-1.11 


+0.32 


-0.67 


-0.83 


+0.22 


+0.15 


-0.05 


+0.17 


-0.59 


-0.74 


+0.24 


-0.68 


-0.77 


+0.30 


-0.37 


-0.52 


+0.20 


-0.16 


-0.33 


+0.10 


+0.23 


+0.02 


+0.08 


+0.62 


+0.39 


+0.19 


+0.20 


-0.01 


+0.20 


+0.58 


+0.34 


+0.21 


-0.65 


-0.77 


+0.14 


-0.35 


-0.50 


+0.19 


-0.36 


-0.49 


+0.17 


-0.23 


-0.32 


+0.08 


+0.69 


+0.46 


+0.14 


-0.78 


-0.97 


+0.31 


+0.24 


+0.09 


+0.08 


-0.75 


-0.86 


+0.21 


+0.26 


+0.11 


+0.11 
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Table 3. Continued. 



JNO. 


KA(J2UUUJ 


Dec(J2000) 




J 


M 


K 


^ eff 




RV 


[M/HJ 




[a/Fe] 




deg. 


deg. 


mag 


mag 


mag 


K 


kms 












(2) 




(3) 


(4) 




(6) 


(') 




(8) 


(9) 


(10) 


(11) 


211 


276.756279 


-33, 


.788357 


13 


.248 


12.575 


12.430 


4404 


-34, 


.348 


-0.54 


-0.70 


+0, 


,19 




276.892143 


-33 


.753590 


13 


.248 


12.670 


1 ^ A On 

12.487 


/I CCI 

4552 


+21. 


.790 


1 A 1 1 

+0.12 


A AC 

-0.05 


+0, 


,10 


'lie 

215 


276.537753 


-33 


.603622 


13 


.256 


12.649 


12.491 


ACCC 

4555 


-119, 


.873 


1 An 
— 1.49 








216 


276.630558 


-33 


.724087 


13 


.258 


12.652 


12.484 


4526 


+38, 


,871 


+0.36 


+0.18 


+0, 


,10 


217 




-33 


.719536 


13 


.258 


12.650 


12.506 


4586 


+53, 


,834 


+0.62 


+0.37 


+0, 


,21 


218 




-33 


.793808 


13 


.263 


12.630 


12.449 


4421 


+25, 


,411 


+0.58 


+0.33 


+0, 


,16 


219 


2/0.491)01 / 


-33 


"7 A 1 AC A 

.791950 


13 


.267 


12.645 


12.524 


4614 


+70, 


,196 


-0.68 


-0.79 


+0, 


,29 


220 


276.786346 


-33, 


C AO i^A 

.598869 


13 


.267 


12.685 


12.496 


4531 


+36, 


,370 


+0.40 


+0.20 


+0, 


,18 


221 


276.488605 


-33, 


.676846 


13 


.277 


12.662 


12.508 


4544 


-68, 


,765 


-0.96 


-0.98 


+0, 


,34 


222 


276.765113 


-33, 


.566185 


13 


.280 


12.582 


12.447 


4371 


+89, 


.637 


-0.84 


-0.94 


+0, 


,11 


ZZi 


276.689571 


-33 


.850536 


13 


.283 


1 1 CA/C 

12.596 


11 AO'i 

12.483 


A A CI 

4452 


+63, 


.721 


1 A 1 A 

+0.19 


1 A nei 

+0.06 


+0. 


,12 


224 


276.530160 


-33 


.772533 


13 


.283 


12.630 


12.546 


4630 


-159, 


.749 


+0.10 


-0.06 


+0. 


,12 


TIC 


276.490013 


-33, 


OCA AAC 

.854095 


13 


.287 


1 1 iCC /I 

12.654 


1 1 CIO 

12.528 


A C£0 

4568 


-2, 


.227 


1 A An 

+0.47 


1 A 10 

+0.28 


+0, 


,10 


226 


276.441487 


-33 


.833378 


13 


.287 


12.601 


12.440 


4344 


+86, 


,733 


-0.55 


-0.71 


+0, 


,21 


227 


Z/D.40U5JU 


-33 


.865948 


13 


.295 


12.693 


12.545 


4593 


-30, 


,132 


+0.42 


+0.20 


+0, 


,12 


228 


276.567606 


-33 


.881664 


13 


.305 


12.729 


12.573 


4641 


+5, 


,293 


+0.17 


-0.05 


+0, 


,12 


229 


0"7*C /lOOTIO 

z/d.486z36 


-33 


1 A y1 O-^ 

.744823 


13 


.308 


12.659 


12.507 


4457 


-93, 


,491 


-0.90 


-0.89 


+0, 


,02 


230 


'^Ty' ^7/^1 rVCT 

276.701957 


-33, 


A r\C'S 

.724052 


13 


.317 


12.670 


12.492 


4390 


-15, 


,167 


+0.44 


+0.21 


+0, 


,17 


231 


276.521504 


-33, 


.882675 


13 


.318 


12.634 


12.517 


4454 


+87, 


.725 


-0.85 


-0.94 


+0, 


,20 


232 


^Ai 1 /TO 

276.701163 


-33 


.812206 


13 


.321 


1 '1 TAI 

12.791 


1 CAO 

12.593 


4650 


+46, 


.717 


A 1 /I 

—0.24 


A OA 

-0.39 


+0, 


,16 


233 


276.742927 


-33 


.655430 


13 


.327 


12.768 


12.577 


A coo 

4588 


-33, 


.817 


■ A /(I 

+0.42 


1 A '^A 

+0.20 


+0, 


,20 


234 


OTiC 1 TOO 1 

276.817331 


-33, 


.599476 


13 


.328 


12.681 


1 1 col 
12.582 


A CAO 

4598 


+3, 


,010 


AIT 

-0.17 


A /( 

-0.34 


+0, 


,06 


235 


Z/D.DOl /9/ 


-33 


TC A /'AO 

.754608 


13 


.331 


12.760 


12.590 


4618 


+2, 


,381 


+0.36 


+0.18 


+0, 


,14 


236 


z/o.ol 1041 


-33 


.747356 


13 


.339 


12.698 


12.552 


4485 


-36, 


,405 


-0.64 


-0.84 


+0, 


,42 


237 


Z/d.d436Z4 


-33 


/' AAO 

.638008 


13 


.349 


12.771 


12.601 


4601 


-128, 


,142 


-1.17 


-1.26 


+0, 


,35 


238 


z/o. V ly 


-33 


.863255 


13 


.363 


12.721 


12.560 


4443 


-48, 


,485 


+0.26 


+0.09 


+0, 


,15 


239 


C CI A AO A 

Z/O.M4464 


-33 


.B939U2 


13 


.367 


12.706 


12.596 


4534 


-6, 


,349 


-0.66 


-0.74 


+0, 


,18 


240 


276.494178 


-33, 


.865170 


13 


.374 


12.671 


12.545 


4384 


-45, 


.410 


-0.45 


-0.61 


+0, 


,22 


^At 

241 


2/6.691228 


-33 


.652210 


13 


.381 


12.740 


12.621 


A CiCO 

4563 


+163, 


.130 


A nc 

-0.75 


-0.82 


+0. 


,21 


242 


OTiC T/TAAOA 

276.769980 


-33 


.622349 


13 


.400 


1 TOC 

12.785 


1 1 con 

12.587 


A A 1A 

4420 


+37, 


.049 


+0.28 


I n A 
+0.14 


+0, 


,08 


243 


276.445539 


-33 


.820004 


13 


.401 


12.802 


12.662 


4626 


-201, 


,762 


-0.51 


-0.68 


+0, 


,34 


245 


Z/D.4^64U(J 


-33 


.680145 


13 


.410 


12.870 


12.674 


4637 


-0, 


,726 










246 


'Sir' '70'^OAC 

276.732205 


-33 


.735016 


13 


.413 


12.722 


12.633 


4506 


-20, 


,965 


-0.26 


-0.40 


+0, 


,19 


247 


z/o.ojz/oU 


-33 


.666142 


13 


.415 


12.859 


12.647 


4542 


-12, 


,734 


+0.63 


+0.34 


+0, 


,22 


248 


z/o.ooo /zz 


-33 


.563961 


13 


.423 


12.782 


12.703 


4680 


-91, 


,605 


-0.82 


-0.94 


+0, 


,39 


249 


276.874535 


-33 


.706005 


13 


.432 


12.805 


12.675 


4564 


+20, 


,002 


+0.22 


+0.05 


+0, 


,05 


TCI 

251 


£i 01 /IT 1 


-33, 


OA 1 1 A 

.891129 


13 


.436 


12.813 


12.703 


4630 


-100, 


.233 


-0.18 


-0.29 


+0, 


,09 


252 


276.590433 


-33, 


.824047 


13 


.447 


1 1 c 

12.815 


12.625 


A 'inn 

4399 


+14, 


.280 


A 1 

-0.16 


A OT 

-0.37 


+0, 


,14 


ICO 

253 


27o.o9oo54 


-33, 


i^t\cn 

.720577 


13 


.450 


1 '1 OAO 

12.802 


12.707 


A£no 

4608 


-12, 


.354 


1 A '^0 

+0.28 


■ A 1 /( 

+0.14 


+0, 


,12 


254 


276.560443 


-33 


TOZn A/I 

.786304 


13 


.457 


1 '1 AAC 

12.905 


1 /CCA 

12.659 


A A iCO 

4463 


+47, 


.473 


■ A OT 

+0.37 


1 A 1 A 

+0.19 


+0, 


,15 


256 


276.885368 


-33, 


A lis /"I 

.847267 


13 


.459 


12.843 


12.729 


4637 


+47, 


.665 


-0.61 


-0.74 


+0, 


,28 


OCT 

257 


Zib.olSZZb 


-33, 


T> AT AA 

.739399 


13 


.464 


12.896 


1 n A A 
12. /44 


A ana 

4676 


+53, 


,227 


1 A C 1 

+0.51 


1 A OA 

+0.30 


+0, 


,21 


258 


276.509126 


-33, 


y' ^ IS Cil 

.638287 


13 


.466 


12.795 


12.648 


4414 


+20, 


,249 


+0.22 


+0.01 


+0, 


,25 


259 


2/D.9U3420 


-33, 


T A y1 /I 1 

.794418 


13 


.475 


12.812 


12.641 


4362 


+0, 


,338 


+0.09 


-0.07 


+0, 


,07 


260 


276.856050 


-33, 


.636570 


13 


.477 


12.791 


12.647 


4375 


-99, 


,944 


-0.31 


-0.39 


+0, 


,09 


261 


276.597652 


-33 


.706051 


13 


.478 


12.889 


12.705 


4530 


-56, 


,387 


+0.25 


+0.06 


+0, 


,14 


262 


2/D.D2/210 


-33, 


.776421 


13 


.481 


12.854 


12.717 


4552 


-31, 


,979 


-0.74 


-0.91 


+0, 


,47 


263 


276.659752 


-33, 


.951805 


13 


.489 


12.884 


12.727 


4553 


-5, 


.578 


+0.04 


—0.12 


+0, 


,09 


264 


'^n£L C/ICO/TA 

276.545369 


-33 


i:\'"MS c\ ✓'A 

.932869 


13 


.493 


1 1 OOA 

12.880 


1 no 1 

12.781 
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4701 


-61, 


.467 


A '1C 

-0.25 


A 00 

-0.38 


+0. 


,17 


266 


276.533760 


-33 


.787659 


13 


.505 


1 '1 on A 

12.874 


1 1 nA A 

12.744 


A CiCI 

4562 


-68, 


.790 


A nc 

-0.75 


A 00 

-0.83 


+0, 


,21 


267 


276.616984 


-33, 


.586750 


13 


.506 


12.954 


12.792 


4699 


-101, 


,202 


+0.59 


+0.35 


+0, 


,18 


268 


2/o.494o /D 


-33 


.647800 


13 


.515 


12.919 


12.794 


4681 


+71, 


,564 


+0.08 


-0.07 


+0, 


,18 


269 


2/0. /o4459 


-33 


C A/" OA 

.589680 


13 


.518 


12.940 


12.784 


4635 


+69, 


,289 


-0.04 


-0.15 


+0, 


,14 


270 


2/d.d1492o 


-33 


AO A'^ 

.938892 


13 


.526 


12.948 


12.808 


4680 


+65, 


,396 


-0.69 


-0.78 


+0, 


,23 


271 


2/d.4o5499 


-33 


y' 1 

.869881 


13 


.526 


12.857 


12.692 


4372 


+29, 


,413 


+0.66 


+0.39 


+0, 


,06 


272 


2/o.o2/432 


-33 


.753338 


13 


.529 


12.833 


12.770 


4559 


-104, 


,732 


-1.17 


-1.37 


+0, 


,38 


273 


276.785331 


-33, 


.659843 


13 


.538 


12.873 


12.710 


4381 


+223, 


.568 


-0.55 


-0.66 


+0, 


,12 


274 


276.842104 


-33 


.732456 


13 


.541 


12.852 


12.714 


4381 


-21, 


.443 


-0.28 


-0.51 


+0, 


,10 


275 


276.765440 


-33 


.581825 


13 


.554 


12.999 


12.817 


4626 


+68, 


.660 


-0.57 


-0.67 


+0. 


,20 


276 


276.728640 


-33, 


.949757 


13 


.557 


12.951 


12.756 


4446 


-121, 


,711 


-0.71 


-0.85 


+0, 


,22 


277 


276.505794 


-33, 


.774723 


13 


.561 


12.941 


12.782 


4515 


-51, 


,159 


+0.33 


+0.05 


+0, 


,16 


278 


276.714377 


-33, 


.724106 


13 


.581 


13.007 


12.860 


4672 


+ 16, 


,094 


+0.29 


+0.09 


+0, 


,12 


279 


276.619384 


-33, 


.610786 


13 


.591 


12.968 


12.842 


4596 


+50, 


,290 


-0.65 


-0.73 


+0, 


,19 


280 


276.699742 


-33, 


.801800 


13 


.600 


13.064 


12.896 


4722 


-22, 


.286 


-0.43 


-0.62 


+0, 


,37 
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276.632739 


-33 


.920406 


13 


.601 


12.931 


12.843 


4565 


-74, 


.670 


+0.19 


-0.02 


+0. 
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Table 3. Continued. 
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RA(J2000) 

deg. 



Dec(J2000) 
deg. 



J 

mag 



H 
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K 
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RV 
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Table 3. Continued. 
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Table 3. Continued. 
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Table 3. Continued. 
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* 



Notes. Meaning of the columns: (1): Target number in this programme; (2): Right ascension from 2MASS; (3): Declination from 2MASS; (4), (5), 
and (6): 2MASS H-, and AT-magnitude; (8): Effective temperature; (7): Heliocentric radial velocity; (9): General metallicity [M/H]; (10): Iron 
abundance [Fe/H]; (11): or-element abundance [a/Fe]; (12): Flag for foreground star candidates. 



